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ABSTRACT 


A  mathematical  model  has  been  developed  to  simulate  the 
interfuel  substitution  process  in  Canada.  This  process  is 
represented  by  the  consumer  when  he  enters  the  energy 
marketplace  to  bargain  for  a  fuel  form  to  satisfy  his  energy 
needs . 

The  model  was  used  to  analyze  the  past  behavior 
patterns  of  energy  consumers  and  to  project  some  possible 
future  scenarios  of  energy  consumption  under  differing 
hypotheses . 
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CHAETEB  I 


INTRODUCTION 

The  Department  cf  Energy,  Mines  and  Resources  of  the 
Canadian  government  has  already  clearly  stated  (4)  the  need 
for  a  comprehensive  understanding  of  the  dynamics  of  energy 
supply  and  demand  in  Canada.  They  correctly  point  out  that 
this  understanding  must  include  knowledge  regarding  the 
interactions  that  the  energy  sector  has  within  a  wide  range 
cf  complex  Canadian  social  and  economic  structures.  From 
this  perspective  the  Department  cf  Energy,  Mines  and 
Resources  has  prcpcsed  a  qualitative  framework  cf  important- 
factors  to  be  considered  when  formulating  energy 
policy  {  4)  1 2 . 

The  necessity  of  clearly  specifying  a  framework  and  the 
corresponding  components  for  analysis  is  well  recognized 
here.  However,  the  following  thesis  is  based  on  the 
conviction  that  an  equally  essential  and  useful  ingredient 
in  the  analysis  methodology  is  the  utilization  cf  'large- 
scale12  mathematical  models,  pa rticular ily  in  instances  such 
as  energy  policy  formulation  where  many  of  the  important 
components  are  easily  quantified.  The  functions  that  'large- 
scale'  mathematical  models  perform  include: 


1  see  figure  1.1 

2  large-scale  models  are  defined  here  as  multiple,  open 
system  models. 
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FIGURE  I. 
DEPARTMENT 


1  -  QUALITATIVE  FRAMEWORK  WITHIN  WHICH  T 

OF  ENERGY,  MINES  AND  RESOURCES  PROPOSES 
FORMULATE  CANADIAN  ENERGY  POLICY. 
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-  Since  it  is  a  relatively  simple  task 
to  divide  large-scale  mathematical 
models  intc  sub-components,  these  models 
facilitate  easier  component  by  component 
analysis  within  the  framework  of  a  fcicaa 
perspective  . 

-  The  quantification  of  variables 
required  by  mathematical  models  produces 
mere  concise  definitions  of  abstract 
mental  concepts. 

-  The  quantitative  results  produced  fcy 
large-scale  mathematical  models  yield  an 
understanding  of  the  relative  importance 
of  interrelated  factors  and  the  impacts 
that  input  changes  may  cause. 

-  In  the  particular  case  of  energy 
policy  analysis,  it  is  vital  that  the 
available  information  concerning  future 
expected  quantities  such  as  energy 
supplies,  demand  and  prices  be  as 
quantitatively  accurate  as  possible. 

Therefore,  the  primary  objectives  of  this  thesis  are  to 
expand  the  information  base  required  to  formulate  energy 
policy  and  to  illustrate  the  use  of  * large-sca le 1 
mathematical  models  as  valuable  tools  in  this  frccess. 


With  these  objectives  in  mind,  a  mathematical  model 
describing  the  interfuel  substitution  process  in  Canada  was 
developed.  Only  the  behavior  of  the  consumer  is  depicted  in 
this  model  wherein  interfuel  substitution  refers  to  the 
process  in  which  consumers  are  engaged  when  they  enter  the 
marketplace  to  bargain  for  energy  forms  to  satisfy  their 
energy  requirements.  This  process  is  a  major  factor  in 
determining  energy  consumption  patterns  and  must  be 
understood  before  policy  decisions  affecting  energy 


consumers  can  be  made. 


E.  ■ 


» 

. 
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CHAPTER  II 


LITERATURE  REVIEW 

Not  until  recently  have  economists  and  government 
p cl icy  analysts  become  concerned  with  problems  related  to 
interfuel  substitution,  This  statement  is  fully  supported  by 
the  fact  that  an  insignificant  quantity  of  literature  on 
this  topic  was  published  before  1970. 

Energy  problems  in  the  first  half  cf  this  century  were 
greatly  simplified  by  the  facts  that  there  was  generally  an 
abundant  supply  of  cheap  energy  and  the  environmental 
impacts  cf  energy  projects  were  relatively  small. 
Technological  changes  since  1900  have  contributed  to  a 
continual  expansion  and  diversification  of  the  Canadian 
energy  base.  In  addition  tc  coal  being  a  valuable  energy 
source,  these  technological  changes  have  resulted  in  crude 
oil,  natural  gas,  and  electricity  also  being  important 
energy  forms  available  to  the  ultimate  users  cf  energy  in 
Canada.  Even  the  producers  of  electricity  have  a  wider 
choice  of  primary  energy  fcrms  available  to  them  as  a  result 
of  these  technc logical  changes.  The  electricity  producers  of 
today  may  choose  from  hydro,  natural  gas,  crude  oil,  coal  or 
nuclear  power  as  primary  sources  of  energy  for  electricity 
gene ration. 

It  was  the  introduction  of  these  alternative  fuel  forms 
that  gave  rise  tc  the  first  real  occurences  of  interfuel 
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competition  in  Canada-  However,  the  situation  that  existed 
until  very  recently  may  be  described  as  a  "consumer's 
market"  where  energy  producers  were  competing  with  each 
ether  for  a  share  of  the  market.  This  mentality  still 
existed  in  1968  (9). 

In  the  late  1S6C's  the  emergence  of  environmental 
concerns  added  a  new  dimension  to  the  already  complex  nature 
of  the  interfuel  ccmjietiticn  process.  Canadians  began  to 
acquire  a  consciousness  cf  the  impacts  that  resource 
development  (including  energy  resources)  and  subsequent 
consumption  have  on  cur  natural  environment.  There  can  be  no 
question  that  environmental  concerns  have  affected  the 
development  and  consumption  patterns  of  energy  resources. 

New  circumstances  since  1S72  have  fundamentally  altered 
the  situation  regarding  energy  resources  in  Canada.  It  is 
new  recognized  that  Canada  is  entering  a  period 
characterized  by  the  scarcity  of  low-cost  energy.  There  is 
an  increasing  pressure  tc  re-examine  our  future  energy 
requirements  and  tc  define  mere  clearly  the  criterion  that 
should  be  employed  in  satisfying  our  energy  demands. 

Included  in  the  decision  criterion  there  must  be  the 
consideration  of  the  technological  restrictions,  comparative 
efficiencies  (19,33)  and  costs,  and  the  environmental 
impacts  of  fuel  utilization.  Manifestations  of  these  new 
concerns  and  issues  are  available  in  the  form  cf  recent 
conferences  and  government  studies  (3,14,15,16,23). 


. 


. 
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Despite  the  fact  that  in  the  last  decade  there  have 
teen  an  increasing  proliferation  of  studies  and  reports 
addressing  the  issues  of  energy  demand  and  interfuel 
substitution,  the  theoretical  foundations  of  these  endeavors 
must  still  be  considered  in  their  infant  stages.  Gonzalez 
(9)  and  Fose  (2  2)  have  analyzed  and  discussed  the 
determinants  and  the  consequences  of  interfuel  competition 
from  the  perspective  of  classical  demand  theory  as  outlined 
by  Samuelson  and  Scott  (24) .  That  is,  through  the 
identification  of  the  price  of  fuels  and  their  substitutes, 
consumer  tastes  and  consumer  incomes,  it  is  possible  to 
qualitatively  trace  the  impacts  that  price,  taste  or  income 
changes  will  have  on  the  consumption  of  particular  fuels, 
their  substitutes,  and  their  cc  if  pi  aments. 

Important  contributions  towards  understanding  the 
utilization  of  energy  resources  have  been  made  by  studies 
applying  an  economic  analysis  to  the  environmental  problems 
(6,14).  However,  these  developments  have  been  cf  a  more 
general  nature  relating  to  the  classification  cf  a  broad 
spectrum  cf  ’resources'  and  the  specification  cf  their 
individual  characteristics. 

In  addition,  substantial  theoretical  advances  have  been 
made  through  the  use  cf  mathematical  simulation  techniques. 
Seari  (2£)  has  described  several  energy  models  which  have 
teen  developed  in  the  United  States.  Foell  (8)  and  Palmedo 
(18)  have  demonstrated  how  the  use  of  the  mathematical 
models  can  guide  the  process  of  energy  research  and  assist 


in  policy  f or mulatic n .  Although  the  Behrens  (2)  model  of 
natural  resource  utilization  does  not  discuss  energy  usage 
per  se,  his  model  structure  does  identify,  in  a  very 
aggregated  way,  the  critical  factors  to  be  examined  in  the 
analysis  of  ary  natural,  ncnrene wable  resource.  Gn  the  other 
hand,  Baughman  (1)  and  Khazzcom  (11,12)  have  loth  conducted 
studies  on  very  specific,  well-defined  aspects  of  energy 
problems.  Martin  Baughman  has  constructed  a  model  of  the 
interrelationship  cf  energy  supply,  demand  and  price  in  the 
United  States.  His  work  in  the  area  cf  interfuel 
substitution  by  consumers  is  one  of  the  first  cuantitative 
approaches  to  this  subject.  Khazzoom  is  currently  in  the 
ficcess  cf  constructing  several  independent  single-equation 
models  of  energy  demand  in  Canada.  These  models  describe 
sector-by-sect cr  demand  for  individual  energy  sources.  It  is 
expected  that  Khazzcom  will  interrelate  these  models  once  he 
has  constructed  equations  for  all  the  sectors  and  all  the 
energy  forms. 

A  very  different  approach  to  energy  demand  problems  has 
been  presented  by  Collins,  Harren  and  Uozniak  (5).  They 
consider  the  possibility  of  inter-sector  cooperation,  rather 
than  competition,  in  the  utilization  of  fuels.  As  a  result, 
they  have  developed  a  scheme  for  the  optimal  allocation  of 
crude  oil,  natural  gas,  and  coal  among  the  power,  plastics 
and  steel  industries. 

Since  the  field  cf  energy  research  within  its  present 
context  is  sc  new  a  literature  review  is  never  complete.. 


. 


* 


* 

Several  times  each  year  new  ideas,  approaches  and  results 


are  emerging. 


This  thesis  was  completed  with  the  understanding  that 
the  above  mentioned  works  constituted  the  major 
contributions  to  energy  demand  research  thus  far. 


.  : 


CHAPTER  III 


GENERAL  DESCRIPTION  AND  DEVELOPMENT  OF  THE  MODEL 

Three  basic  stages  are  involved  in  developing  a  ’large- 
scale'  mathematical  model.  In  the  first  stage  the 
theoretical  foundations  of  the  model  are  outlined  in  a 
qualitative  manner.  It  is  in  this  stage  of  model  development 
that  the  model  assumptions,  levels  of  aggregation,  and  the 
major  variables  and  their  interrelationships  are  identified. 
The  second  stage  of  the  model  development  consists  of 
quantifying  the  variable  interrelationships  in  a  suitable 
mathematical  form.  Finally,  in  the  third  stage,  the  validity 
cf  the  mathematical  model  is  tested  and  the  model  as  it 
exists  is  either  accepted,  accepted  subject  tc  improvements, 
or  rejected.  Although  the  stages  of  model  development  have 
boon  described  as  occurring  in  a  single  sequential  process, 
these  stages  are  generally  followed  in  a  repetitive  fashion 
as  the  model  builder  converges  on  a  valid  model;  i.e. 
theory,  mathematical  representation,  validity,  theory, 
mathematical  representation,  etc.. 

The  development  of  the  Canadian  Interfuel  Substitution 
Model  will  be  described  with  reference  to  the  three  stages 
cf  model  development. 
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JU  THEORY 

The  basic  concepts  required  to  understand  the  economics 
cf  energy  demand  and  interiuel  substitution  have  already 
been  outlined  and  expanded  upon  by  Khazzcom  (11)  and 
Eaughman  (1).  The  development  cf  these  theories,  as 
extensions  cf  the  classical  demand  theory,  originated  around 
1963  (4)  .  In  cider  tc  save  the  reader  time,  a  summarized 
account  cf  the  theory  will  he  presented  here,. 

The  total  demand  for  energy  resources  exists  in  two 
fundamental  forms  -  captive  demand  and  free  demand. 

Captive  demand  is  that  part  of  the  total  demand  which 
is  immobilized  by  past  commitments.  These  past  commitments 
may  include  investment  in  single-fuel  transportation  systems 
(i.e.  natural  gas  pipeline)  or  single-fuel  conversion 
equipment  (i.e.  electric  stove).  Because  of  these 
commitments,  the  consumer  is  net  in  a  position  of 
ccn side ring  interfuel  substitution. 

On  the  other  hand,  free  demand  is  that  portion  of  the 
total  demand  which  is  free  from  any  commitments  and  may  be 
satisfied  by  any  number  of  energy  fuel  forms  depending  on 
the  economic  and  technological  conditions  prevailing  at  the 
time.  The  free  demand  of  any  given  sector  is  equal  to  the 
sum  of  the  replacement  demand  and  incremental  demand  for 
that  particular  sector.  Beplace.ment  demand  originates  from 
consumers  who  were,  in  the  past,  committed  to  a  particular 
fuel.  However,  these  consumers  are  no  longer  tied  to  any  one 
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fue-l  and  they  are  new  entering  the  energy  marketplace  in 
search  of  possible  substitute  fuels.  The  consumers  who 
possess  replacement  demand  behave  as  if  they  were  in  the 
energy  marketplace  for  the  first  time.  That  is,  they 
consider  all  fuels  as  possible  substitutes  including  the 
particular  fuel  that  they  utilized  prior  to  entering  the 
marketplace.  Incremental  demand  is  that  portion  cf  total 
demand  made  up  cf  new  consumer  needs  or  incremental  growth 
in  total  demand.  Because  cf  the  sensitivity  of  free  demand 
tc  market  conditions,  it  is  often  referred  to  as  the  ’market 
sensitive’  demand. 

The  ratio  cf  free  demand  tc  captive  demand  for  energy 
resources  in  Canada  has  generally  been  very  small  in  the 
past.  This  situaticn  is  partly  explained  by  the  fact  that 
the  costs  cf  conversion  equipment  have  been  relatively  high 
compared  to  the  fuel  costs.  As  a  result,  consumers  have  been 
comimitted  by  large  investments  to  their  particular 
conversion  equipment  and  they  are  unable  or  unwilling  to 
make  spontaneous,  short-range  decisions  concerning  fuel 
choices.  However,  it  is  possible  that  the  ratio  cf  free 
demand  tc  captive  demand  will  increase  in  the  future  if  fuel 
prices  escalate  at  a  faster  rate  than  conversion  equipment 
costs . 


The  criteria  that  energy  consumers  with  free  demand 
employ  in  committing  their  demand  to  a  particular  fuel  form 
are  identical  tc  those  identified  in  classical  demand 
theory.  These  criteria  consist  of  a  relative  comparison  of 
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fuel  prices  and  the  prices  cf  their  complements  (i.e.  fuel- 
consuming  appliances)  and  a  consideration  of  the  non¬ 
economic  requirements  (i„e.  tastes)  of  the  consumer.  A 
contemporary  example  of  the  effects  that  changing  tastes  may 
have  on  the  interfuel  substitution  process  is  the  rise  of 
environmental  concerns  producing  demands  for  cleaner  energy 
fcrtns.  In  addition,  it  must  be  remembered  that  consumers  may 
diminish  thsii  free  demand,  and  hence  their  total  demand, 
for  energy  fuels  by  substituting  labour  or  capital  for  their 
requirements.  Examples  cf  this  would  be  the  substitution  of 
higher  quality  insulation  for  space  heating  cr  the 
substitution  cf  a  bicycle  for  an  automobile. 

This  concludes  the  summary  cf  the  basic  theoretical 
concepts  required  to  analyze  the  demand  for  energy  and  the 
interfuel  substitution  process.  As  the  results  and 
conclusions  are  discussed  in  later  chapter,  mere  detailed 
theoretical  observations  will  be  made. 
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E.  1  HE  MCEEL 

A  mathematical  model  ^hich  simulates  the  in ter fuel 
substitution  process  for  energy  demand  in  Canada  has  been 
developed.  The  model,  with  its  levels  of  aggregation, 
framework,  and  mathematical  structure  has  been  adopted  in 
its  basic  form  frciu  the  work  dene  by  Baughman  fl).  The  only 
major  structural  changes  that  have  been  made  tc  the  Baughman 
mcdel  relate  to  the  nature  of  the  '  distribution 
multipliers'.  The  distribution  multipliers  are  a  set  of 
mathematical  relationships  which  simulate  consumer  behavior 
in  the  energy  marketplace.  Ihese  multipliers  determine  what 
proportions  cf  the  total  free  demand  for  energy  will  be 
satisfied  by  each  cf  the  fuel  forms  available.  A  more 
detailed  description  of  the  distribution  multipliers  will  be 
given  later  1 . 

Energy  demand  in  the  model  is  assumed  to  arise  from 
feur  sources  -  residential  and  commercial  space  heating, 
industrial  process  heating,  transportation,  and  electricity 
generation.  In  turn,  these  energy  demands  are  assumed  to  be 
satisfied  by  seme  combination  of  the  four  "end-use"  energy 
forms  available  -  crude  oil,  ccal,  natural  gas,  and 
electricity.  Crude  oil,  ccal  and  natural  gas  are  referred  to 
as  primary  energy  forms  since  they  are  converted  directly 
from  their  natural  occurring  forms  into  useful  energy.  On 
the  other  hand,  electricity  is  considered  a  secondary  energy 


1  refer  to  Appendix  E 
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form  because  it  is  produced  from  primary  energy  sources  such 
as  crude  oil,  coal,  natural  gas,  uranium  and  hydro. 

For  convenience,  all  energy  demands  and  their 
corresponding  supplies  are  treated  in  the  model  in  their  BID 
(British  Thermal  Unit)  equivalent  forms.  The  conversion 
factors  are  presented  in  Appendix  A. 

Figure  III.1  graphically  depicts  the  model  with  its 
boundaries,  levels  of  aggregation,  and  interrelationships  of 
variables. 

The  demand  sectors  represented  in  the  model,  with  the 
exception  of  the  electricity  generating  sector,  all  have 
identical  mathematical  structures.  The  electricity 
generating  sector  will  later  be  treated  separately.  Thus, 
the  primary  consuming  sectors  {transportation,  residential 
and  commercial,  and  industrial)  each  have  within  them  three 
irajcr  su  t-ccmponent  s  : 

1.  determination  of  sector  free  demand 

2.  calculation  of  distribution  multipliers 

3.  determination  of  current  sector  demand  for  each 

fuel. 


In  the  process  cf  simulating  the  sectors*  interfuel 
substitution  behavicr  the  model  calculates,  in  a  sequential 
fashion,  the  values  cf  the  three  sub-component s  for  each 
tine  period  and  returns  the  new  values  in  preparation  for 
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the  calculations  at  the  next  tiis  period. 

Since  tbe  mathematical  notation  for  subscripts,  which 
is  employed  in  the  following  explanations,  is  slightly 
irregular  it  will  new  be  briefly  expounded  upon.  The 
variable  subscripts  are  represented  by  m  and  n  and  they 
identify  two  basic  characteristics  of  the  variables  with 
which  they  are  associated.  The  subscript  m  identifies  the 
demand  sector  (residential-commercial ,  transportation, 
industrial  or  electricity  generating)  which  the  variable 
refers  to  while  the  subscript  n  identifies  the  particular 
energy  source  (coal,  oil ,  gas  or  electricity)  which  the 
variable  refers  to. 

1.  EETEEHOATION  Of  SECTCE  EEEE  DEMAND 

As  outlined  in  section  III. A,  the  sector  free  demand  is 
composed  of  fcwc  components  -  replacement  demand  and 
incremental  demand.  According  tc  its  respective  components, 
sector  free  demand  (market  sensitive  demand)  may  be 
described  as  fellows: 

ES  Cm  (t )  =  I Dm  (t)  +  HDm  (t)  (III.1) 

I  Era  (t)  =  EGm  x  TDm  (t)  (III. 2) 

EEm(t)  =  [ raCLPIn  (t)  x  mDn  (t)  ]  (III.  3) 

where  m  -  the  demand  sector  in  question 

n  =  the  energy  fuel  in  question 

flSDj  =  the  market  sensitive  demand  fer  sector  ra  at 
time  t. 

I Dm  =  the  incremental  demand  for  sector  m 
EDm  -  the  replacement  demand  for  sector  m 
EGm  =  the  rate  of  growth  of  energy  demand  for 
sector  m 

TDm  =  the  total  energy  demand  of  sector  m 
mCIEIn  -  the  commitment  liberation  rate  for 


. 
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consumers  of  fuel  n  In  sector  m 
mDn  =  the  demand  for  fuel  n  by  sector  m 

That  is,  the  market  sensitive  demand  of  sector  m  (MDSj) 
equals  the  sum  cf  the  incremental  demand  (IDj)  and  the 
replacement  demand  (EDm)  for  that  sector* 

The  incremental  demand  for  the  total  sector  (IDa)  is 
simply  the  sector  rate  cf  growth  (RGm)  multiplied  by  the 
total  sector  demand  {TDj)  in  the  given  year.  Fcr  each 
primary  consuming  sector  the  rate  cf  growth  (RGm)  is  an 
exogenous  input  and  equal  to  the  average  rate  cf  growth  that 
existed  for  that  sector  in  the  period  1945-70.  Therefore, 
the  incremental  demand  for  each  primary  sector  is  an 
exogenous  input  to  the  model. 

The  replacement  demand  for  sector  m  (HD  a)  is  determined 
by  the  summation  cf  the  products  of  the  consumer  commitment 
liberaticn  rates  fcr  fuel  n  (mCLBTn)  and  the  current 
consumption  cf  fuel  n  in  sector  m  (jDn) .  The  concept  of 
commitment  liberation  rates  as  employed  in  the  'Canadian 
Interfuel  Subst it u ticn  Model*  is  very  broad.  Since  the 
notion  of  the  commitment  liberation  rate,  as  used  here,  is 
synonymous  with  the  rate  at  which  consumers  unlock  past 
commitments  to  enter  the  marketplace,  it  is  obvious  that 
possible  reasons  for  unlocking  past  commitments  would 
include  physical  depreciation  of  equipment,  economic 
benefits,  technological  obsolescence  and  convenience.  In  the 
model,  the  consumer  commitment  liberation  rates  (mCLRIn)  are 


exogenous  inputs. 
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2.  CALCULATION  CP  DISTRIBUTION  MULTIPLIERS 


The  function  that  the  distribution  multipliers  perform 
in  each  demand  sector  is  to  determine  what  fraction  of  the 
sector's  free  demand  is  satisfied  by  each  fuel  form  in  each 
time  period. 

There  are  several  functional  forms  which  may  be  used  to 
represent  the  distribution  multipliers.  Regardless  of  which 
functional  form  is  accepted,  there  are  two  general 
constraints  which  must  be  satisfied  by  the  distribution 
multipliers  within  each  demand  sector.  First,  all  of  the 
distribution  multipliers  must  be  non-negative  and  secondly, 
the  sum  cf  all  the  distribution  multipliers  within  any  given 
demand  sector  must  egual  1.0  for  every  time  period. 

Only  two  functional  forms  for  the  distribution 
multipliers  will  be  discussed  here.  These  twc  forms  will  be 
referred  to  as  the  Baughman  functional  form  and  the 
Mcdif ied-Khaszccm  functional  form.  The  attempts  to  derive 
distribution  multipliers  for  the  Canadian  Interfuel 
Substitution  Model  assumed  that  these  functions  were 
strictly  economic  functions  in  that  they  could  adequately  be 
represented  by  price  variables  alone. 

The  Baughman  functional  form1,  which  was  discarded  as 


1  Martin  Baughman,  Dynamic  Energy  System  Modeling  - 
Interfuel  Competition,  Report  #72-1,  Energy  Analysis  and 
Flaming  Group,  School  of  Engineering,  M.I.T.,  Cambridge, 
Mass.  August  1972,  p.  63. 
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an  acceptable  ferns  fer  the  Canadian  Interfuel  Substitution 

Model,  is  justified  primarily  cn  the  basis  of  one  of  its 

characteristics.  The  distribution  multipliers  were 

constrained  by  definition  to  be  non-negative  since 

Baughman’s  functions  were  represented  by  log-linear 

relationships.  The  structure  of  the  Baughman  distribution 

multipliers  for  each  demand  sector  are  as  follows: 

In  (d)  =  (A  x  F)  (III. 4) 

where  In  (d)  -  natural  logarithm  of  d 

d  =  n  x  1  matrix  of  distribution  multipliers 
n  =  number  of  fuel  forms  available 

=  n  x  n  matrix  of  price  coefficients,  all  values 
are  assumed  tc  be  constant. 

E  =  n  x  1  matrix  of  fuel  prices. 

Be-arranging  equation  III. 4  yields, 

d  =  cEXP  ( Ax  P  )  (III. 5) 

where  c  =  n  x  1  matrix  of  constants. 

EXF(AxE)  =  exponent  of  (AxP) 

In  order  tc  ensure  that  the  sum  of  the  distribution 
multipliers  is  1.0,  Baughman  divided  the  initial  estimates 
by  their  sum  for  each  time  period.  This,  in  effect, 
normalizes  the  initial  estimate  and  guarantees  that  the  sum 
cf  the  distribution  multipliers  will  be  1.0. 

Unfortunately,  the  validity  of  the  Baughman  functional 
form,  as  it  applies  to  the  Canadain  situation,  came  under 
serious  questioning  after  parameter  estimation  studies  were 
conducted.  The  signs  cf  most  of  the  estimated  parameters 
were  different  from  those  that  were  predicted  by  simple 
economic  theory.  In  addition,  many  of  the  parameter 
estimates  were  found  to  be  insignificantly  different  from 
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ze  re  in  a  statistical  sense .  For  these  reasons  it  was 
decided  tc  reject  the  Baughman  functional  form  and  search 
for  a  more  appropriate  fcrm  to  suir  the  Canadian  Interfuel 
Substitution  Model. 

The  work  completed  by  Khazzoom  (11,12)  suggests  that 
the  interfuel  decision  process  is  more  accurately  depicted 
ly  a  comparison  cf  the  relative  fuel  prices  rather  than  a 
comparison  of  the  absolute  fuel  prices  as  suggested  by 
Baughman.  As  a  result  of  Khazzoon's  findings  an  attempt  was 
made  to  represent  the  distibution  multipliers  as  polynomial 
functions  of  the  relative  fuel  prices.  This  functional  form 
is  referred  to  as  the  Modi f ied-Khazzoora  functional  form.  It 
was  clear  from  the  results  cf  the  initial  parameter 
estimates  of  the  Mod if ied-Khazzccm  functional  fcrm  that  the 
decision  criterion  fcr  allocating  free  demand  was  dependent 
cn  more  than  simply  the  relative  prices  of  fuels.  In  spite 
of  this  acknowledgement  it  was  decided  that  a  complete 
specification  cf  the  distribution  multipliers  would  be  a 
major  undertaking  which  cculd  net  be  justified  within  the 
framework  of  this  thesis.  Therefore,  attempts  were  made  to 
represent  the  non- price  variations  in  the  distribution 
iiu ltipliers  with  time  series  variables.  This  approach  proved 
successful  in  most  cases. 

Thus,  the  Modif ied-Khazzcom  distribution  Multipliers 
fcr  the  Canadian  Interfuel  Substitution  Model  possess  the 
following  general  structure: 
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=  aO  +  aj  [PnSATIO(t)  ]  +  a2  [  PnRATIO  (t )  ]2 

+  a3  [TIME]  (III. 6) 

where  mdn  (t)  =  the  distribution  multiplier  for  fuel  n  in 

sector  m 

ak  =  estimated  parameters  of  independent 
variables,  k  •=  0,1,2, 3. 

PnRATIO(t)  =  the  ratic  of  the  price  cf  fuel  n  to 
the  geometric  mean  of  other  fuel  prices. 

TIME  =  time  series  variable. 

As  a  result  cf  the  success  achieved  in  estimating  the 
parameters  for  the  Mcdif ied-Khazzoo®  functional  form  with 
the  additional  time  series  variables,  these  distribution 
multipliers1  were  accepted  for  use  in  the  Canadian  Interfuel 
Substitution  Model.  Cne  unavoidable  drawback  cf  these 
functions  was  the  necessity  cf  imposing,  in  an  ’artificial' 
manner,  the  ncn-negative  and  unity  constraints  cf  the 
distribution  multipliers.  The  non-negative  constraint  was 
satisfied  by  eguating  the  value  of  the  distribution 
multiplier  to  either  its  initial  estimate  or  zero  - 
whichever  is  larger.  The  unity  constraint  was  satisfied  by 
employing  the  same  normalizing  process  which  Eaughman  used. 

3.  DE1ES  MI N ATIC  K  Cl  CUBEENI  SECTOR  DEMAND  EY  EACH  FUEL 

In  the  simulation  cf  the  primary  energy  demand  sectors 
the  determination  of  the  current  sector  demand  for  each  fuel 
represents  the  final  step  for  any  given  time  period.  Here, 
the  cumulative  results  cf  earlier  calculations  are 
summarized  as  the  'updated*  estimates  cf  the  sector-by- 


1  refer  to  Appendix  E  fcr  complete  derivation  cf  the 
Mcdef ied-Khazzccm  functional  form  and  the  corresponding 
statistical  results  cf  parameter  estimation. 
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sector  demands  for  each  energy  form  are  presented. 

The  structural  form  from  which  these  new  estimates  are 
calculated  for  the  primary  energy  demand  sector  is  as 
fellows: 
d[jDn  (t)  ] 

- Z - =  -  [mfiDn(t)  ]  +  [man  ( t )  x  MSDm(t)  ]  (III.  7} 

dt 

where  mDn  =  the  demand  for  fuel  n  by  sector  m 

mEDN  =  the  replacement  demand  for  fuel  n  in 
sector  m 

mdn  =  the  distribution  multiplier  for  fuel  n  in 
sector  nj 

MSDm  -  the  market  sensitive  demand  in  sector  m 

Therefore,  the  instantaneous  rate  of  change  of  demand 
for  fuel  n  by  sector  m  at  time  t  is  equal  to  the  replacement 
demand  for  fuel  n  in  sector  i  subtracted  from  the  proportion 
cf  market  sensitive  demand  which  is  allocated  to  fuel  n  in 
sector  m.  Clarification  cf  this  equation  will  be  given  later 
in  this  chapter  when  a  detailed  description  of  the 
Residential  and  Commercial  demand  sector  equations  will  be 
presented.  By  integrating  this  equation  over  time  it  is 
possible  to  obtain  the  level  of  demand  for  fuel  n  in  each 
sector  m. 

Thus,  the  operation  cf  the  model  as  it  sinulates  the 
interfuel  substitution  process  of  the  primary  demand  sectors 
may  be  summarized  as  follows.  First,  the  values  of  the 
exogenous  varibles  for  each  sector  are  determined,.  The 
exogenous  variables  are  the  total  sector  energy  demand,  rate 
cf  growth  cf  total  sector  demand,  individual  fuel  prices  and 


. 
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the  consumer's  commitment  liberation  rate  for  each  fuel. 

With  this  infer ma tic n  the  model  calculates  the  total  market 
sensitive  jfree)  demand  fer  each  sector;  converts  this  free 
demand  through  the  distribution  multipliers  into  demand  for 
particular  fuel  forms  and  then  estimates  the  current  sector 
demand  for  each  fuel. 

However,  as  mentioned  earlier,  the  electricity 
generating  sector  is  a  special  case  since  it  consumes 
energy,  produces  energy  and  has  more  fuel  forms  available  to 
it  than  any  ether  sector.  The  electricity  generating  sector 
receives,  as  cutput  frem  the  primary  consuming  sectors,  the 
total  and  incremental  demands  for  electricity.  The 
cuantities  of  electricity  produced  from  hydro  and  nuclear 
power  are  given  as  exogenous  inputs  to  the  model.  With  this 
information  the  model  is  able  to  determine  what  proportion 
of  the  free  demand  is  satisfied  by  fossil  fuels.  Then,  the 
calculations  proceed  in  much  the  same  manner  as  for  the 
primary  consuming  sectors  to  determine  the  electricity 
generating  sector's  demand  fer  each  of  the  fossil  fuels. 

Having  completed  a  general  description  cf  the 
structures  of  the  primary  demand  sub-models  and  the 
electricity  generating  sub-mcdel,  a  detailed  description  of 
one  of  the  primary  demand  sub-models  and  the  electricity 
generating  sub -models  will  new  be  given. 

For  illustrative  purposes,  the  residential  and 
commercial  space  heating  sub-model  will  be  listed  and 
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discussed  in  full  detail1.  The  sub-model  listirg  given  below 
will  be  in  mathematical  notation, 

RESIDENTIAL  AND  COEBERCIAL  DEMAND  SOE-MODEL 


The  residential  and  commercial  demand  is  modeled  as  an 
exponential  growth  process.  The  rate  of  change  of  the  sector 
demand  is  given  by 
d[ BCD  (t)  ] 

- —  =  £CDEG  x  BCD  (t)  (III. 8) 

dt 

FCE  (1945)  =  ,38  (III. 9) 

FCDEG  =  .071  (III. 10) 

where  BCD  =  the  level  of  the  residential-c cm mere  la  1 

demand  initialized  at  .38  Q*s2  for  the  year 
19  45. 

RCDGB  =  rate  of  growth  of  8CD  modeled  at  7.14  per 
year  in  the  base  case. 

The  rate  cf  change  cf  the  consumption  of  fuel  n  in  the 
resident ial-ccmmercial  sector  is  given  by 
d£  RCDn  (t )  ] 

- =  - RCDnP.D  (t)  +  [RCDDn(t)  x  RCMSD(t)  j  (III.  11) 

dt 

fer  n  =  W,X,  Y,Z  (coal,  gas,  oil  and  electricity 
respectively) 

where  RCDn  (t)  =  the  level  of  residential-commercial 

consumption  cf  fuel  n  at  time  t. 

RCCnSD(t)  =  residential-commercial  replacement 

demand  for  fuel  n  at  time  t  which  also  equals 
the  rate  cf  decline  of  consumption  of  fuel  n 
if  none  of  the  market  sensitive  demand  is 
supplied  by  fuel  n  at  time  t, 

RCMSD  (t)  =  residential-commercial  market  sensitive 

demand  at  time  t. 

RCEDn  (t)  =  the  residential-commercial  demand 

distribution  factor  multiplying  the  market 
sensitive  demand  for  each  fuel  at  time  t. 


1  see  Appendix  C  for  complete  listing  of  the  me del. 

2  one  Q  corresponds  to  1015  ETU's. 
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The  initial  level  of  consumption  for  each  of  these 
fuels  is 


PCD W  (194  5)  = 

.27  Q*  s 
initi al 

se  ctor 

demand 

for 

coal. 

(III.  12) 

R  C  D  X  (19  4  5)  = 

.02  Q*s 
initial 

sector 

de mand 

for 

natural  gas. 

(III.  13) 

RCDY(1945)  = 

.03  Q  f  s 
initial 

sector 

demand 

for 

crude  oil. 

(III.  14) 

RCDZ  (1945)  = 

.06  Q 5  s 

i  n  i.  t  i  a  1 

sector 

demand 

for 

electricity. 

(III. 1 5) 

The  replacement 

demands 

;  are  given 

by 

BCDnRD  (t)  =  RCDn  (t)  x  RCnB  (III.  16) 

for  n  =  {f,  X,  Yf  Z. 

where  RCDnRD  =  the  resident ial-ccmmercial  replacement 

demand  for  fuel  n , 

RCnR  =  tne  commitment  liberation  rate  for 
consumers  of  fuel  n . 

The  residential-commercial  market  sensitive  demand  is 
the  sum  of  the  incremental  demand  and  replacement  demands; 
it  can  be  written  as 

EC  MS  D  (t )  =  [  RCDRG  x  RCD  (t)  ]  *  ^RCDnRD(t)  (III.  17) 

n 

where  RCMSD(t)  =  residential-commercial  market  sensitive 

demand  at  time  t. 

The  distribution  multipliers  are  calculated  in  a  two- 
step  procedure.  First  the  initial  estimates  are  computed, 
then  they  are  normalized  so  their  sum  is  equal  to  unity. 

This  is  in  line  with  the  assumption  that  total  consumption 
in  each  consuming  sector  is  inelastic.  First,  the  estimates 
are  given  by 
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RCDDnl  (t)  =  anO  -  [  arpl  x  PnRATIO  (t)  ]  +  [ an 2  x  PnRA TIQ  ( t ) 2  ] 

+  [ an3  x  TIME  ]  (III. 1 8) 

for  n  =  W,X,Y,Z. 

where  RCDDnl  =  initial  estimate  of  the  distribution 

multiplier  for  fuel  n  at  time  t, 
an_k  -  estimated  coefficients  for  distribution 
multiplier  of  fuel  n  for  k  =  0,  1,2,3. 

PnRATIO  (t)  =  ratio  of  the  price  of  fuel  n  to  the 

geometric  mean  of  the  fuel  prices  of  the 
other  three  fuels  at  time  t, 

TIME  =  time  in  years  Initialized  at  1945. 

The  sum  of  these  Initial  estimates  is  given  by 
RTOTAL  (t)  =  ^SCDDnl  (t)  '  (III.  19) 

n 

and  the  normalized  distribution  multipliers  are  given  by 
RCDDn(t)  =  RCDDnl  (t)  /  RTOTAL(t)  (III.  20) 

for  n  =  W,X,Y,Z, 


where  RCDDn  (t)  -  the  normalized  distribution  multiplier 

which  is  used  to  allocate  the  free  (market 
sensitive)  demand  at  time  t. 

Therefore,  the  above  are  the  basic  equations  for  the 

residential  and  commercial  demand  sub-model. 

ELECTRICITY  GENERATING  SUB-MODEL 


The  electricity  generating  sub-model  is  quite  similar, 
in  structure,  to  that  just  described  for  the  primary 
consuming  sectors.  However,  besides  generating  electricity 
from  the  fossil  fuels,  the  model  must  take  into  account  the 
role  of  the  nuclear  and  hydro  generation.  The  additional 
electrical  power  generated  by  nuclear  and  hydro  sources  in 
each  time  period  is  an  exogenous  input.  Therefore,  the 
quantity  of  market  sensitive  demand  which  is  satisfied  by 
hydro  and  nuclear  sources  in  each  time  period  is 
predetermined.  Once  the  market  sensitive  demand  is 
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determined  endogenously,  the  quantity  of  electricity  which 
must  be  generated  from  fossil  fuels  is  calculated  by 
subtracting  the  incremental  hydro  and  nuclear  generation 
from  the  total  market  sensitive  demand.  Then,  the  choice  of 
fossil  fuel  mixes  is  made  in  the  same  manner  as  in  the 
primary  demand  sectors.  The  market  sensitive  demand  for  the 
electricity  generating  sector  is  simply  the  sum  of  the 
incremental  demand  for  electricity  of  the  primary  demand 
sectors  and  the  replacement  demand  of  the  fossil  generated 
electricity. 


The  incremental  electricity  demand  can  be  written  as 

ZDG1 (t)  =  ^{increments  in  electrical  consumption  from  each 
primary  consuming  sector)  (III, 21) 

where  ZDG1  (t)  =  the  incremental  electricity  demand  at 

timet. 

The  fossil  fuel  replacement  demands  are  given  by 
ZFnRD(t)  =  ZFn(t)  x  ZFn  B  (III.  22) 

for  n  =  W , X ,  and  Y. 

where  ZFnED  (t)  =  the  replacement  demand  in  electricity 

for  fossil  fuel  n  at  time  t. 

ZFn  (t)  =  the  quantity  of  total  electrical  output 

supplied  by  fuel  n  at  time  t. 

ZFnB  =  the  commitment  liberation  rate  for 

consumers  of  fuel  n  in  the  electricity 
generating  sector. 

The  fraction  of  electrical  output  supplied  by  nuclear 
power  is  assumed  to  be  the  same  as  the  fraction  of  total 
capacity  made  up  of  nuclear  power. 

ZFN  (t)  =  [1.0  -  FCF(t)  ]  x  Z0(t)  (III.  23) 

where  ZFN(t)  =  the  electical  output  produced  from 

nuclear  generation  (Q* s/year)  at  time  t. 

FCF(t)  =  the  fraction  of  capacity  made  up  of 
fossil  fuel  and  hydro  plants  at  time  t. 
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ZO  (t)  -  total  electrical  output  at  time  t. 

The  electrical  output  produced  from  hydro  generation 
(ZFH )  is  exogenously  fed  into  the  model  via  a  Dynamo  table 
function. 


The  fossil  incremental  demand  is  derived  by  subtracting 

the  increments  in  nuclear  and  hydro  output  from  the  total 

growth  in  electrical  output.  This  is  written  as 

ZFFG(t)  =  ZDG1  (t)  -  (AZFN  ( t)  /At)  -  (AZFH  (t )  /At)  (III.  24) 

where  ZFFG  (t)  =  the  growth  in  electrical  output  to  be 

supplied  by  fossil  fuels  at  time  t, 

ZDG1 (t)  =  the  total  growth  in  electrical  output  at 
time  t. 

AZFN/At  -  the  growth  in  output  supplied  by  nuclear 
generation. 

AZFH/At  -  the  growth  in  output  supplied  by  hydro 
generation. 

The  fossil  market  sensitive  demand  is  then  the  sum  of 
the  fossil  incremental  demand  and  the  fossil  replacement 
demand,  or 

ZFSD(t)  =  ZFFG  (t)  +  SZFnBD(t)  (III.  25) 

n 

where  ZFSD  (t)  =  the  fossil  market  sensitive  demand  at 

time  t. 

ZFFG  (t)  -  the  fossil  incremental  demand  at  time  t. 
ZFnRD  (t)  -  the  replacement  demand  of  fuel  n  at 

time  t. 

It  is  possible  that  a  shift  from  fossil  generation  to 
hydro  or  nuclear  generation  may  occur,  in  which  case  ZFFG  (t) 
in  Equation  111.24  would  be  negative. 

The  dynamics  of  the  fossil  fuel  demands  in  electricity 
are  now  given  by  the  same  equations  as  those  for  the  primary 
demand  sectors. 


The  consumption  of  the  primary  fuels  (in  Q’s)  in 
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electricity  generation  is  easily  obtained  front  the  equations 
nlOZ(t)  =  [  ZFn  ( t )  x  BBF  <t)  ]/ (.  0094)  (III. 26) 

for  n  =  v  ,x ,Y. 

where  nTCZ  =  the  consumption  of  fuel  n  in  the 

electricity  sector  at  time  t. 

ZFn  (t)  -  the  output  of  electricity  produced  by 
fuel  n  at  time  t „ 

HBF  { t )  =  the  fossil  heat  rate  in  millions  of  BTU's 

per  Kilowatt-hour  at  time  t.  The  same  heat 
rates  that  Baughman  used  were  assumed  here. 

.0094  =  the  lossless  conversion  rate  in  millions 
of  BTO's  per  kilowatt-hour. 

This  concludes  the  description  of  the  electricity 
generating  sector.  Although  only  the  residential-commercial 
demand  sector  and  the  electricity  generating  sector  were 
described  in  detail  here,  a  complete  listing  cf  the  total 
model  is  available  in  Appendix  C 


The  task  cf  transforming  the  qualitative  theory  of 
interfuel  subs t ituticn  into  a  corresponding  quantitative 
form  is  complete.  The  final  stage  in  the  development  of  the 
Canadian  Interfuel  Substitution  Model  consists  cf  testing 
the  validity  cf  the  mathematical  model. 


i 


the  listing  in  appendix  C  is  given  in  Dynamo  (21)  format. 
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C.  J5CpjI_VALir&TION 

I  he  problem  of  validating  mathematical  models  in  the 
social  sciences  is  very  complicated  and  it  has  received 
little  attention  in  tbs  past® 

First,  we  must  have  a  clear  understanding  cf  the  term 
•validity*.  Many  social  scientists  have  implicitly  accepted 
the  concept  of  validity  that  equates  ’valid’  with  ’true’. 
«hen  applied  to  mathematical  models,  this  approach  has 
several  drawbacks.  It  is  improper  to  consider  these  models 
in  such  strong  terms  as  being  either  true  or  false.  A  social 
model  can  never  be  proven  unconditionally  true.  In  addition, 
this  approach  fails  tc  recognize  the  benefits  gained  from 
constructing  the  model  by  ever-emphasizing  the  final 
results. 

Probably  the  most  satisfactory  definition  cf  the  term 
validity  in  the  context  of  the  Canadian  Inter fuel 
Substitution  Model  is: 

A  model  is  valid  if  it  is  useful  for  a  clearly 
stated  purpose.  And,  for  causal  descriptive 
models,  a  purpose  calls  for  a  model  which 
describes  reality  with  sufficient  accuracy  to  be 
useful  for  investigating  alternative  social 
policies  1 . 

kith  this  definition  in  mind,  an  outline  of  the  model 
validation  procedure  can  be  described.  Generally  speaking. 


1  Eeter  M .  Senge  "Seme  Issues  in  Evaluating  the  Validity  of 
Social  System  Models",  Proceedings  of  the  Summer  Computer 
Simulation  Conference  (Montreal  1973),  p.  1176. 
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two  broad  categories  cf  ircdel  testing  techniques  exist:  the 
tests  cf  whole  ircdel  behavior  (Macro  tests)  and  tests 
related  to  the  structural  components  (Micro  tests).  Within 
each  category  there  are  several  tests  which  may  be  employed. 
An  exhaustive  validation  procedure  was  not  implemented  due 
tc  time  constraints  although  the  tests  conducted  here  are 
considered  sufficient. 


MACFC  TESTS 


i.  Easic  Assumptions  cf  Model  Structure 

Under  this  heading  the  basic  assumptions  underlying  the 
model  structure  are  examined  and  if  they  are  unacceptable 
the  necessary  changes  are  made. 

ASSUMPTION  I:  Energy  consumer  decisions  can  be  adequately 
modeled  when  the  levels  cf  aggregation  are  such  that  there 
exist  four  national  sources  cf  energy  demand 
(transportation,  residential  and  commercial,  industrial  and 
electricity  generation)  which  are  satisfied  by  four  forms  of 
energy  fuels  (crude  oil ,  natural  gas,  coal,  and 
electricity) . 

It  is  obvious  from  the  historical  data  that  Canadian 
energy  consumption  patterns  have  varied  widely  from  one 
geographic  region  tc  another.  However,  the  only  consequence 
cf  this  fact  is  tc  place  limitations  on  the  corclusions  and 
recommendations  that  may  be  derived  from  the  model. 


The  levels  of  aggregation  used  in  the  model  present  no 
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difficulties  for  model  validation  if  the  infer mat ion 
extracted  from  the  ircdel  is  cn  a  corresponding  level  of 
aggregation. 

ASSUMPTION  II:  In  the  model  there  has  been  no  explicit 
consideration  of  the  possibilities  of  technological  changes 
which  may  introduce  new  fuel  terms  (i.e.  solar  energy)  or 
perhaps  even  create  new  demand  sectors.  Although  these 
unaccounted  possibilities  are  part icular ily  important  for 
long  range  projections,  they  are  not  considered  to  be 
significant  for  the  15  year  projection  span  used  in  this 
study. 


MICFC  TESTS 

i.  £ata  Ease  Validation 

The  data  base  for  the  model  was  double  checked  in 
search  of  possible  errors  and  problems  of  inaccuracy  or 
misinterpretations. 

Although  several  problems  dc  exist  in  the  raw  data,  the 
necessary  corrections  and  adjustments  were  made  to  make  the 
data  and  the  model  structure  compatible.  A  full  description 
of  this  data  base  and  the  associated  problems  is  given  in 


Appendix  A. 
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ii.  Basic  Assumptions  of  Structural  Components 
ASSUMPTION  I:  The  commitment  liberation1  rates  cf  consumers 
are  exogenous  inputs  into  the  model  and  they  have  assumed 
values  since  nc  data  are  available  for  these  values.  With 
the  exception  cf  two  cases  (IHWE/TRWB)  ,  these  values  are 
considered  constants  and  they  have  rather  small 
values  (<.20)  according  to  intuitive  judgement. 

Generally,  it  is  agreed  that  the  longer  a  fuel  has  been 
consumed  in  significant  quantities,  the  higher  will  be  the 
potential  commitment  liberation  rate  of  consumers  using  this 
fuel.  For  this  reason,  the  commitment  liberation  rates  of 
coal  consumers  have  been  very  high  while  those  of  gas 
consumers  are  still  relatively  lew. 

The  commitment  liberation  rates  used  in  the  model  are 
believed  to  adeguately  represent  past  behavior.  However, 
problems  do  arise  when  providing  future  expected  values  of 
these  commitment  liberation  rates  for  model  projections. 
ASSUMPTION  II:  The  fuel  prices  at  the  production  sites  {coal 
mines,  gas  fields,  etc.)  are  sufficiently  correlated  with 
the  retail  fuel  prices  sc  that  they  may  be  used  as  the 
prices  that  consumers  respond  to. 

Of  all  the  assumptions  made,  this  one  appears  to  be  the 
most  g ues tic nable .  large  fluctuations  in  retail  fuel  prices 


1  as  defined  earlier,  the  commitment  liberation  rates  of 
consumers  refer  to  the  rates  at  which  consumers  unlock  their 
past  committments  to  particular  fuel  forms. 
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dc  exist  fro®  one  geographic  region  to  another.  In  addition, 
there  are  several  cases  where  it  is  a  foreign  fuel  price 
that  consumers  ate  responding  tc  {USA  coal  imported  to 
Ontario,  foreign  crude  oil  imported  to  Eastern  Canada)  . 
Although  this  later  point  is  important  when  examining  the 
determinants  of  price  levels  and  price  changes,  it  is  much 
less  significant  when  examining  only  the  impact  of  these 
prices  on  demand  since  the  foreign  fuel  and  domestic  fuel 
prices  and  fluctuations  have  been  very  similar  since  1945. 

In  spite  cf  the  undesirability  of  this  assumption 
concerning  fuel  prices,  it  was  made  in  order  tc  simplify  the 
model.  As  long  as  this  simplification  does  not  lead  to  gross 
inaccuracies  or  misuses  cf  the  model  it  is  considered 
acceptable. 

AS  SUMPTION  III:  The  nuclear  and  fossil  heat  ccrversion  races 
are  the  same  as  these  used  by  Baughman  (1) .  Although  they 
may  be  different  for  Canada,  the  difference  is  assumed  to  be 
negligible. 

iii.  'Gccdness  cf  Fit'  Tests 

The  area  cf  model  validation  which  has  received  the 
most  attention  in  the  past  is  that  associated  with  'goodness 
of  fit*  tests.  It  is  under  this  heading  that  mcdel  builders 
describe,  generally  in  a  quantitative  manner,  the  accuracy 
with  which  their  mcdel  is  able  tc  reproduce  past  data. 

Several  quantitative  tests  that  perform  this  task  have 
teen  developed.  The  'Canadian  Int-arfuel  Substitution  Model' 
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was  subjected  to  only  three  of  these  tests. 


MEAN  PERCENTAGE  ERROR  ~  This  measure  is  simply  the  mean  of 
all  the  percent  errors  that  exist  between  the  actual  data 
and  the  model  predictions.  A  problem  with  this  measure 
exists  in  that  a  declining  function  will  invariably  have 
large  percentage  errors  at  its  lower  values  even  if  the 
absolute  errors  remain  unchanged.  This  produces  an  over¬ 
weighting  of  the  errors  at  lower  values  of  a  function.  The 
mathematical  form  of  this  measure  is  : 


1 

^MODELm  -  ACTUAL^ 

MPE  = 

-  *1 

' 

}.T 

ii 

pi  ACTUALm  ) 

(III. 27) 


where  N  =  the  number  of  data  sample  points. 


EQUATION  VARIANCE  -  The  equation  variance  is  another  measure 
of  the  dispersion  of  the  actual  data  about  the  estimated 
model  equation.  The  mathematical  form  of  this  measure  is: 


VAR 
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N 


MODELm 


ACTUALm) 2 


where  N  =  the  number  of  data  sample  points. 


(III. 28) 


There  is  a  tendency,  with  this  statistic,  to  give  more 
weighting  to  the  larger  model  values  since  the  absolute 
error  is  likely  to  be  larger  there. 


R-SQUARED  -  Although  this  test  statistic  is  really 
associated  with  the  Ordinary  Least  Squates  technique  of 
estimating  parameters  for  linear  equations,  it  can  be 
applied  outside  of  this  scope  with  the  appropriate 
limitations.  The  purpose  of  this  measure  is  to  indicate  what 
proportion  of  the  variation  of  the  data  about  its  mean  is 
described  by  the  estimated  equation.  A  value  of  R -Squared  = 
1.0  would  indicate  a  perfect  model  fit  while  deviations  from 
this  value  would  indicate  less  than  perfect  fits.  The 
mathematical  form  of  this  measure  is: 

N 


(MODELm  -  ACTUALm) 2 
(ACTUALm  -  MEANACT)2 


(III. 29) 


where  N  =  the  number  of  data  sample  points, 

MEANACT  =  mean  of  the  actual  data  for  the  period  1  945- 


70. 


The  R-Squared  measure  for  goodness  of  fit  is  less  suitable 
for  rapidly  changing,  nonlinear  (i.e.  exponential)  functions 
since  that  test  statistic  under  these  conditions  will  tend 
to  give  an  over-optimistic  impression  of  the  goodness  of 
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Since  it  is  recognized  that  each  of  the  test  statistics 
described  above  has  its  own  advantages  and  disadvantages,  it 
is  believed  that  an  accurate  measure  of  the  gccdness  of  fit 
would  be  obtained  by  using  all  three  of  these  statistics. 

The  results  of  the  gccdness  of  fit  tests  are  presented 
separately  for  the  smoothing  functions  which  were  used  to 
simulate  the  exogenous  data  of  the  model  and  the  endogenous 
variables  of  the  model.  It  is  important  to  make  this 
distinction  when  one  is  analyzing  the  source  of  model 
errors.  In  the  case  cf  the  smoothed  exogenous  data  of  the 
model,  the  gccdness  cf  fit  tests  are  applied  tc  those 
functions  which  appear  as  estimated  time  series  functions. 
The  estimation  errors  that  exist  in  these  functions  have 
direct  effects  on  the  final  results  of  the  model 
predictions.  Cn  the  ether  hand,  the  errors  that  arise  in  the 
endogenously  determined  variables  are  the  results  cf  errors 
in  the  smoothing  functions  of  the  exogenous  data  coupled 
with  the  inability  cf  the  model  to  fully  explain  consumer 
behavior . 

Table  I2I.1  and  Figures  III. 2  -  III. 3  present  the 
results  cf  the  goodness  cf  fit  tests  for  the  major  smoothing 
functions  of  the  exogenous  data  of  the  model. 

Table  III. 2  and  Figures  III. 4  -  III. 15  present  the 
results  cf  the  gccdness  cf  fit  tests  for  the  major 
endogenous  variables  of  the  model,. 


. 


. 


. 
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At  this  point  it  must  be  mentioned  that  an  explicit 
optimization  procedure  was  not  used  to  obtain  the  best  fit 
of  the  model  to  the  data.  Rather,  a  trial  and  error  approach 
was  employed  to  minimize  the  errors  and  slight  improvements 
are  likely  possible. 

Figures  III. 16  -  III. 20  present  the  output  of  the  Ease 
Case  which  was  net  subjected  to  goodness  of  fit  tests.  The 
term  'Ease  Case'  refers  to  the  results  of  the  model 
simulation  for  the  period  1945-7C. 


1ABLF  III.1 

'Goodness  of  Fit'  Statistics  for  the  Smoothing  Functions  of 

the  Exogenous  Data  in  the  Model 


Ee  pendent 

Equation 

Equation 

Eg  ration 

Variable 

R2 

Variance 

Percentage 

RCD 

.9937 

.001  9 

2.67 

I  HD 

.9854 

.0029 

3.03 

TED 

.9150 

.0024 

4.64 

. 


, 


. 


« 


. 

• 

PRIMARY  DEMAND  SECTOR  (Q) 


2  *  5C004- — 


2  f  coco— 


1  *■  5000- 


V  - 

+  - 

A  - 

x  - 


Total.  Residential-Commercial 
Demand  (Actual  Data) 

Total  Resident ial-C  ommerci a 1 
Demand  (Model  Prediction) 

Total  Transportation  Demand  (Actua 1 
Data) 

Total  Transportation  Demand  (Model 
Predict  ion) 


v 

+ 


V 


V 

+ 


^  4- 


V 


4* 


+ 


V 

+ 


FIGURE  III. 2  -  TOTAL  ENERGY  DEMAND  BY  RESIDENTIAL-COMMERCIAL 
AND  TRANSPORTATION  SECTORS,  BASE  CASE  1945-7C. 
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FIGURE  III. 3 


TOTAL  ENERGY  DEMAND  BY  INDUSTRIAL  SECTOR , 
BASE  CASE  1945-70. 


IAEIE  III. 2 


Goodness  of  Fit.*  Statistics  for  the  Endogenous  Variables  in 

the  Model 


Dependent 

Equation 

Variable 

F2 

TPEC 

.  5936 

TED 

.  9  955 

ECCW 

.9513 

RCDX 

.9957 

EC  DY 

.  5649 

ECDZ 

.  9940 

IHE» 

.  54  68 

IHBX 

.  9881 

IHDY 

.  5456 

IHDZ 

.  9821 

TEEW 

.9610 

TED  Y 

.  9548 

ZFN 

.  8742 

ZFH 

.  8501 

ZFX 

.  64  34 

ZFY 

.73  10 

WD 

.9719 

XD 

.  9921 

YD 

.9816 

ZD 

.9930 

Equation 

Equation 

Variance 

Percentage 

.0099 

2.  84 

.0065 

1.73 

.0003 

14.  06 

.0000 

7.  34 

.0024 

8.13 

.0003 

4.  85 

.0002 

4.83 

.0003 

25.  31 

.0010 

11.40 

.0007 

2.77 

.0006 

141.02 

.003  9 

6-  06 

.0000 

----  i 

.0010 

43.  79 

.0001 

- i 

.0001 

44.76 

.0026 

3.74 

.0007 

13.  95 

.0089 

5.  48 

.0013 

2.30 

Mean 

Error 


1  Meaningful  estimates  are  not  available  since  the 
percentage  error  is  indeterminant  for  a  large  number  of  the 
sample  data  points. 
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FIGURE  III. 5  -  TOTAL  OF  HAND  FOR  NATURAL  GAS  (XD),  BASE  CASE 

1 945-70. 
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FIGURE  III. 7  - 


TOTAL  CONSUMPTION  OF  PRIMARY  ENERGY  FORMS 
(TP EC)  ,  BASE  CASE  1  94  5-70. 
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III. 8  -  RESIDENTIAL-COMMERCIAL  ENERGY  DEMAND,  BAGS 

CASE  1945-70. 
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FIGURE  III. 9  -  RESIDENTIAL-COMMERCIAL  ENERGY  DEMAND,  BASE 

CASE  1945-70. 
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FIGURE  III. 1C  -  INDUSTRIAL  ENERGY  DEMAND,  BASE  CASE  1945-70. 
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FIGURE  III. 11  -  INDUSTRIAL  ENERGY  DEMAND,  BASE  CASE  1945-7C. 


TRANSPORTATION  DEMAND  (Q 


H 


sH 


H 

VI 


H 


' — I 


H 


E  (YEARS) 


FIGURE  III. 12  -  TRANSPORTATION  ENERGY  DEMAND,  BASE  CASE 

1945-70. 
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FIGURE  III .13- 


TOTAL  DEMAND  FOR  ELECTRICITY  (ZO)  ,  BASE  CASE 
1945-70. 
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FIGURE  III. 14  -  ELECTRICITY  GENERATION  FRO M  NUCLEAR  AND 
NATURAL  GAS,  BASE  CASE  1945-70. 
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FIGURE  III. 15  -  ELECTRICITY  GENERATION  FROM  COAL  AND  CRUDE 

OIL,  BASE  CASE  1945-70. 
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FIGURE  III. 16  -  FUEL  PRICES,  BASE  CASE  1945-70. 
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FIGURE  III. 17  -  RELATIVE  PRICES  OF  FUELS,  BASE  CASE  1945-70. 
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FIGURE  Til. 18 


TOTAL  MARKET  SHAPES  OF  PRIMARY  FUELS, 
CASE  1945-70. 
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FIGURE  III. 19  -  MARKET  SHARES  OF  PRIMARY  FUELS  IN 
ELECTRICITY  GENERATING  SECTOR,  BASE  CASE  1945-70. 
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FIGURE  III. 20 


QUANTITY  OF  ELECTRICITY  PRODUCED  FROM  HYDRO, 
BASE  CASE  1945-70. 
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iv,  Per t urbatio n  Tests 

A  useful  exercise  in  validating  mathematical  models  is 
to  subject  tie  model  tc  large-scale  shocks  and  to  observe 
the  results  for  anticipated  cr  unexpected  behavior.  If  the 
results  of  the  perturbation  tests  are  unexplainable  and 
obviously  incorrect,  the  validity  of  the  model  may  then  be 
questioned . 


Perturbation  tests  sere  conducted  on  the  Canadian 
Interfuel  Substitution  Model  by  introducing  radical  changes 
in  fuel  prices,  fossil  capacity  fraction  and  the  quantity  of 
electricity  produced  from  hydro.  Although  the  results  are 
net  presented  here,  the  mcdel  passed  the  perturbation  tests 
since  nc  unacceptable  results  were  observed. 

v.  Eealism  cf  Mcdel  Parameters 

Are  the  values  assigned  tc  parameters  within  the  mod^l 
realistic  in  their  eccncmic  interpretation? 

Certainly  all  those  parameters  associated  with  the 
distribution  multipliers  are  realistic  since  they  were 
statistically  estimated  from  actual  data.  The  estimation 
procedure  included  examination  for  the  expected  signs  of  the 
parameters  where  important  in  order  to  determine  whether  or 
net  they  correspond  to  economic  theory. 

As  mentioned  earlier,  the  commitment  liberation  rates 
correspond  tc  assumed  values  which  are  already  intuitively 
correct  in  their  magnitude. 


. 


. 
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It  is  new  believed  that  the  model  has  successfully  been 
demonstrated  tc  be  valid  according  to  the  definition  of 
validity  used  here.  Tc  deny  that  inadequacies  exist  in  the 
mcdel  would  be  unrealistic.  The  model  must  be  understood 
with  all  of  its  strengths  and  weaknessess  and  any  statements 
or  conclusions  must  be  made  in  the  light  of  the  capabilities 
of.  the  model. 


, 

' 

. 

CHAPTER  IV 


ANALYSIS  AND  RESULTS 


In  order  to  exercise  and  assess  t 
'Canadian  Interfuel  Substitution  Model 
examination  of  the  model  will  attempt 
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A.  GENERAL  OBSERVATIONS 


The  primary  demand  sector  growth  rates  have  been  larger 
in  Canada  than  the  United  States  for  the  period  1945-70. 
There  are  two  major  factors  contributing  to  Canada's  larger 
growth  rate  in  energy  consumption: 

1. )  The  Canadian  average  annual  population  growth 
rate  has  been  significantly  larger  than  that  for 
the  USA.  It  is  generally  accepted  that  the  total 
energy  consumption  of  a  country  is  closely 
correlated  with  its  population  size. 

2. )  The  level  of  Canadian  energy  consumption  per 
capita  has  been  lower  than  that  in  the  USA  and 
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growing  at  a  faster  rate.  Although  a  controversial 
theory,  it  has  been  hypothesized  (32)  that  the 
energy  consumption  per  capita  in  a  growing  economy 
follows  an  'S'  curve  over  time,  eventually 
reaching  a  saturation  level.  Proponents  of  this 
theory  point  to  the  fact  that  the  American  energy 
consumption  per  capita  is  growing  at  a  declining 
rate  and  they  believe  that  it  is  nearing 
saturation. 


Table  XV*  1  contains  the  comparative  data  cf  Canadian 
and  American  energy  consumption  patterns. 


TABLE  IV.  1 

Growth  Rates  of  Energy  Consumption  for  Canada  and  USA 


(19  45-70) 


Population 

Primary  Demand 

Annual 

Growth  Rate 

Growth 

Sec  to  r 

Average  Per 

cf  Annual 

Rate 

Growth  Sates 

Capita  Energy 

Average  Per 

(%/Year) 

(%/Year) 

SCD  I HD  TRD 

Consul pt ion 
(BTUxlO6) 

Capita  Energy 
Consumption 
(%/Year) 

Canada  2.26 

7.16  4.55  3.05 

181.2 

2.72 

USA  1  . 74 

3.74  2.47  2.70 

217.7 

1.  31 

Within  the  primary  demand  sectors  of  Canada  there  have 
existed,  since  1945,  relatively  smooth  and  definite  trends 
with  regards  to  the  preferences  for  particular 
energy  forms1. 


Coal  has  made  a  transition  from  being  the  most  demanded 


1  refer  to  Figures  IV. 3  and  IV. 4,  the  time-varying  fuel 
market  shares. 
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FIGURE  IV. 1 


ELECTRICITY  GENERATION  FROM  PRIMARY  FUELS , 
BASE  CASE  1945-70. 
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FIGURE  IV. 2  -  ELECTRICITY  GENERATED  FROM  HYDRO,  BASF  CASE 

1945-70. 
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FIGURE  IV. 3  -  PRIMARY  FUEL  MARKET  SHARES  IN  ELECTRICITY 
GENERATING  SECTOR,  BASE  CASE  1945-70. 
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FIGURE  IV. 4  -  PRIMARY  FUEL  MARKET  SHARES  OF  TOTAL  CANADIAN 
ENERGY  CONSUMPTION,  BASE  CASE  1945-^0  . 
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fuel  form  in  1945  to  being  the  least  demanded  fuel  form  in 
197C.  It  is  estimated1  that  in  1945  coal  constituted  50%  of 
all  the  primary  energy  forms2  utilized  in  Canada.  By  197C 
the  market  share  of  coal  had  declined  to  less  that  10%. 
Several  factors  have  contributed  to  this  transition:  the 
prices  of  competitive  fuels  declined  relative  to  the  price 
of  coal;  serious  transportation  and  handling  problems 
existed  with  coal;  and  environmental  concerns  discriminated 
against  ccal  as  a  'dirty1  fuel. 

During  the  period  1945-55,  the  importance  of  crude  oil 
as  a  primary  source  of  energy  grew  very  rapidly.  However, 
since  1955  the  more  rapid  growth  in  natural  gas  consumption 
has  prevented  crude  oil  from  increasing  its  market  share. 

The  electricity  generating  sector  has  exhibited 
sporadic  behavior  in  the  past  with  regards  to  its  demands 
for  particular  fuel  forms3.  This  sector  is  a  public  utility 
in  most  provinces  of  Canada.  The  monopolistic  operation  of 
the  electric  utilitites  has  enabled  this  sector  to  operate 
or.  a  very  large  scale  with  adherence  to  long  range 
objectives  and  planning.  As  a  result,  the  following 
characteristics  may  be  attributed  to  the  electricity 
generating  sector: 

1.)  Based  on  environmental,  technical,  and  cost 
considerations,  hydro  and  nuclear  generated 


1  refer  to  'Base-Case'  market  shares. 

2  the  primary  fuels  are  coal,  natural  gas,  crude  oil, 
and  nuclear  (W,X,Y,H,  and  N  respectively). 

3  see  Figures  IV. 1  to  IV. 3,  electricity  market  shares. 


hydro. 
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electricity  can  be  justified  only  when  operating 
at  full  or  near-full  capacity.  Therefore,  large- 
scale  hydro  and  nuclear  projects  are  brought  to 
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The  above  characteristics  explain  a  large  measure  of 
the  fluctuating  behavior  of  the  electricity  generating 
sector.  In  spite  of  its  sporadic  and  seemingly  random 
behavior,  the  electricity  generating  sector  is  gradually 
reducing  its  hydro  market  share1.  This  is  necessary  in  view 
of  the  fact  than  demand  for  electricity  is  steadily  growing 
while  the  potential  for  hydro  power  is  nearing  its  upper 
limit.  Although  the  lost  market  share  of  hydro  power  has 
been  captured  by  fossil  power  in  the  last  decade,  it  is 
expected  that  nuclear  power  will  become  a  major  source  of 
electrical  energy  in  the  near  future. 


The  preceding 
factors  contribute 
preferences.  These 
equipment  costs,  a 
convenience  of  use 
roles  and  relative 


discussion  recognizes  that,  several 
to  the  determination  of  a  consumer's  fuel 
factors  include  fuel  prices,  conversion 
vailability,  technological  constraints  and 
.  Unfortunately,  the  quantification  of  the 
impacts  that  all  these  factors  have  in 


1  see  Figure  IV. 3, 


electricity  market  shares 
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the  total  decision-making  process  is  extremely  difficult.  In 
fact,  a  meaningful  discussion  of  the  non-price  factors  would 
require  a  regional  analysis  of  consumer  behavior. 

E.  DEMAND  ELASTICITIES 

There  are  primarily  two  characteristics  of  an  energy 
demand  sector  which  determine  the  magnitude  of  that  sector’s 
demand  elasticity  for  any  given  fuel.  These  two 
characteristics  are: 

1. )  the  ratio  of  free  (market  sensitive)  demand 
which  is  satisfied  by  fuel  a  to  the  total  demand 
for  that  fuel. 

2. )  the  elasticity  of  the  free  demand  which  is 
satisfied  by  fuel  n  with  respect  to  the  relative 
price  for  that  fuel. 

In  the  event  that  the  ratio  of  free  demand  which  is 
satisfied  by  fuel  n  to  the  total  demand  for  this  fuel  is 
very  small  tor  any  given  sector,  it  is  virtually  guaranteed 
that  the  price  elasticity  of  total  demand  for  this  fuel  will 
also  be  small.  This  is  true  since  the  base  demand  for  any 
fuel  is  by  definition  perfectly  inelastic.  However,  if  the 
ratio  of  free  demand  to  total  demand  for  a  fuel  is  large, 
then  the  price  elasticity  of  the  total  demand  for  the  given 
fuel  in  this  sector  is  dependent  on  the  elasticity  of  the 
free  demand  which  is  satisfied  by  the  particular  fuel  in 
question . 

Therefore,  in  analyzing  the  price  elasticity  of  demand 
for  each  fuel  in  eaca  sector,  the  following  data  will  be 
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presented  for  each  consumer  category  for  the  period  1945-70: 

-  ratio  of  free  demand  which  is  satisfied  by  fuel  n  to 
total  demand  for  that  fuel. 

-  price  elasticity  of  free  demand  which  is  satisfied  by 
each  fuel  form. 

-  price  elasticity  of  total  demand  for  each  fuel  form. 


Before  presenting  the  results,  the  derivation  of  the 
elasticity  functions  will  be  discussed.  The  market  sensitive 
demand  which  is  satisfied  by  any  one  fuel  is  represented  in 
the  model  by  the  following  function: 

mMSDn(t)  =  raDDn(t)  x  raMSD(t)  {IV.  1) 

for  m  =  BCD, I  HD,  and  TED* 
n  =  w ,X, Y,  and  Z* 

if  a  =  2 F  {electricity  generating  sector) 
then  Equation  IV, 1  is  true  for  n  =  W,X,  and  Y. 


where 


mMSDn(t)  •=  the  quantity  of  free  demand  which  is 
satisfied  by  fuel  n  in  sector  m  at  time  t* 
mDDn  (t)  =  the  distribution  multiplier  for  fuel  n 

in  sector  m  at  time  t. 

jhMSD  (t)  =  the  total  free  demand  in  sector  m  at 
time  t. 


From  equation  IV, 1  it  is  possible  to  derive  a  general 
expression  for  the  elasticity  of  the  free  demand  which  is 
satisfied  by  any  given  fuel* 

As  a  result  of  the  manner  in  which  the  price  variables 
enter  into  the  model,  the  price  elasticities  may  be 
considered  with  respect  to  two  different  price  functions  -- 
elasticity  with  respect  to  absolute  price  and  elasticity 
with  respect  to  relative  price  (price  ratio) .  Preliminary 
studies  of  the  elasticities  with  respect  to  absolute  prices 
discouraged  further  work  with  this  approach  since  the 


elasticity  functions  were  very  complex  and  difficult  to 
manipulate*  In  contrast  to  this  approach,  the  elasticities 
with  respect  to  the  relative  fuel  prices  are  rather  simple 
and  manageable  functions.  Therefore,  elasticity  as  it  will 
be  used  here,  will  new  be  defined  as  the  sensitivity  of 
consumer  demand  for  a  particular  fuel  to  changes  in  the 
relative  price  of  that  fuel. 

The  elasticity  of  free  demand  which  is  satisfied  by 
fuel  n  is  thus  given  by: 

d[ mMSDn  { t)  ]  PnRATIO  ( t) 

E[  mMSDn  (t)  ]  =  - = - x - 

d[  PnRATIO  (t)  ]  mMSDn  (t) 

d[  mDDn  { t )  x  mMSD(t)]  PnEATIO(t) 

= - x -  (IV ,  2) 

d[ PnRATIO  (t)  ]  mMSDn  (t) 

where  E[ mMSDn  (t)  ]  =  the  elasticity  of  the  free  demand  in 

sector  a  which  is  satisfied  by  fuel  n  (mMSDn (t) ) 
with  respect  to  the  relative  price  of 
fuel  n  (PnRATIO  (t)  )  at  time  t. 

Since  the  total  market  sensitive  demand  is  not 
dependent  on  fuel  price  in  the  model,  equation  IV. 2  reduces 
to : 


mMSD  (t)  x  d[mDDn(t)  ]  PnRATIO  (t) 

E[mMSDn(t)]  = - x -  (IV. 3) 

d[  PnRATIO  (t)  ]  mMSDn  (t) 


One  example  of  the  general  structure  of  these 
elasticity  functions  may  be  obtained  by  solving  the 
derivative  of  the  distribution  multiplier  with  respect  to 
the  relative  fuel  price.  This  has  been  done  for  the 
industrial  distribution  multiplier  for  electricity.  The 
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resulting  elasticity  function  is: 


E[  IHMSDZ  (t)  ] 


-0.0112  X  [  ITOTAL  (t)  -IHDD21  (t)  ]  x  1HMSD  (t) 


where 


ITOTAL  (t)  2 
PZEATIO  (t) 


x 


[  IHDDZ  (t)  x  IHHSD  (t)  ] 


(17.4) 


IHMSDZ (t)  =  the  quant 

industrial  sector  whi 
electricity  at  time  t 
ITOTAL (t)  =  the  sum  o 

all  the  distribution 
industrial  sector. 
IHCDZ1  (t)  =  the  initi 

distribution  multipli 
IHDDZ  (t)  =  the  normal 
for  electricity  in  t h 
IHHSD  (t)  =  the  total 
sector  at  time  t. 
-0.0112  =  the  coeffic 
electricity  in  the  in 
electricity  distribut 
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ch  is  satisfied  by 

» 

f  the  initial  estimates  of 
multipliers  for  the 

al  estimate  of  the  industrial 

er  for  electricity. 

ized  distribution  multiplier 

e  industrial  sector. 

free  demand  in  the  industrial 

lent  of  the  relative  price  of 
itial  equation  for  the 
ion  multiplier. 


Slight  variations  of  this  structure  exist  for  other 
elasticities  depending  on  the  particular  form  of  the 
distribution  multiplier  involved. 

Two  basic  characteristics  of  the  elasr ici ties  of  the 
distribution  multipliers  are  clear  from  the  equation: 

1.)  the  constant  term  in  the  elasticity  function 
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elasticities  must  take  into  account,  the  stochastic 
nature  of  this  constant  term, 

2.)  the  elasticity  of  the  distribution  multiplier 


1  see  Appendix  B  for  the  estimation  results. 
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increases  in  absolute  value  as  the  fuel  market 
share  of  free  demand  for  that  sector  declines. 

This  result  is  justified  on  the  basis  that  as  one 
particular  fuel  captures  increasingly  larger 
proportions  of  the  market/  there  is  a  real  decline 
in  the  number  and  extent  of  ‘effective’  competitor 
fuels.  However/  this  feature  of  the  elasticities 
can  only  be  considered  from  the  standpoint  of 
short-run  responses. 


It  is  known  by  definition  that  the  base  demand  for  any 
given  fuel  is  perfectly  inelastic.  Therefore/  we  may  easily 
obtain  the  elasticity  of  total  demand  for  fuel  n  in  sector  m 
by  modifying  the  elasticity  of  free  demand  satisfied  by  fuel 
n  in  sector  m.  This  is  done  by  multiplying  the  elasticity  of 
free  demand  by  the  ratio  of  free  demand  satisfied  by  fuel  n 
to  the  total  demand  for  fuel  n. 


d[  mMSDn  (t)  J  Pn  RATIO  (t )  mMSDn  (t) 

£[mdn(t)  ]  = - - - x - x -  (IV.  5) 

d[  PnFATIO  (!)  ]  mMSDn(t)  ffidn(t) 

mMSD(t)  x  d[  mDDn  (1)  ]  PnSATIO  (t) 

= - x -  (IV, 6) 

d[ Pn RATIO  (t)  ]  mdn  (t) 


Thus,  the  two  elasticity  functions  in  which  we  are 
interested  in  are  represented  mathematically  by  equations 
IV . 3  and  IV. 6,  Equation  IV. 3  determines  the  elasticity  of 
the  free  demand  which  is  satisfied  by  fuel  n  in  sector  m. 
Equation  IV. 6  determines  the  elasticity  of  the  rotal  demand 
for  fuel  n  In  sector  m. 


These  elasticity  functions  will  now  be  discussed  in 
regards  to  the  demand  sectors  of  the  Canadian  Interfuel 


Substitution  Model, 
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EES IDE  NT IAL  AND  COMMERCIAL  SECTOR! 

Coal  demand  in  this  sector  has  been,  since  1945,  much 
more  sensitive  to  price  changes  than  rhe  demand  for  any 
ether  fuel.  This  fact,  coupled  with  the  declining 
residential-commercial  coal  marker  share,  results  in  the 
elasticity  of  free  demand  which  is  satisfied  by  coal  being 
the  largest  of  all  the  demand  elasticities  in  this  sector. 
The  elasticity  cf  the  free  demand  which  is  satisfied  by  coal 
continues  to  increase  as  the  coal  market  share  of  free 
demand  declines. 

As  the  coal  consumers  in  the  residential-commercial 
sector  begin  to  substitute  other  fuels  for  coal  there  is  a 
marked  decline  in  the  ratio  of  free  demand  which  is 
satisfied  by  coal  to  the  total  demand  for  coal.  This  trend 
influences  the  elasticity  of  total  demand  for  coal  by 
greatly  reducing  it  as  compared  to  the  elasticity  of  free 
demand  for  coal.  Despite  this  fact,  the  trends  of  the 
elasticities  of  total  and  free  demand  for  coal  are 
remarkably  similar. 

Crude  oil  demand  in  the  residential-commercial  sector 
has  been  completely  insensitive  to  price  changes  since  1945. 
Thus,  the  price  elasticity  of  the  free  demand  which  is 
satisfied  by  crude  oil  has  been  zero.  Although  the  ratio  of 
free  demand  to  total  demand  for  crude  oil  has  been  quite 
dynamic  since  1945,  the  elasticity  of  the  total  demand  for 

see  Figures  IV. 5  to  IV. 


i 
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this  fuel  is  nonetheless  zero. 

The  elasticity  of  the  free  demand  committed  to 
electricity  has  been  steadily  declining  since  1945  as  the 
proportion  of  free  demand  satisfied  by  electricity 
increases.  The  elasticity  of  total  demand  for  electricity  in 
this  sector  has  generally  been  low  (<.05)  since  the  ratio  of 
free  demand  to  total  demand,  for  this  fuel  has  been  small. 

The  behavior  of  the  elasticities  of  the  demand  for 
natural  gas  have  been  very  similar  to  that  of  electricity. 
INDUSTRIAL  SECTOR1 

In  the  past,  the  coal  distribution  multiplier  for  the 
industrial  sector  has  been  independent  of  price  and  hence 
the  demand  elasticities  have  been  zero. 

All  other  fuels  in  this  sector  have  had  relatively 
inelastic  fuel  demands  and,  with  the  exception  of  natural 
gas,  these  elasticities  have  been  quite  stable  since  1945. 

The  comparatively  large  elasticity  of  total  demand  for 
natural  gas  is  explained  in  part  by  the  large  ratio  of  free 
demand  to  total  demand  that  these  consumers  possessed.  This 
is  particularly  true  for  the  period  1945-55. 

It  may  be  of  interest  to  note  at  this  time  the 
relationship  between  the  ratio  of  free  demand  to  total 
demand  for  fuel  n  and  the  length  of  time  in  which  this  fuel 


1  see  Figures  IV, 9  to  IV. 11. 
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has  been  consumed  in  the  given  sector, 
being  consumed  for  the  first  time  by  a 
result  in  a  very  large  ratio  of  free  d 
The  magnitude  of  the  demand  ratio  will 
which  the  given  sector  adopts  this  fue 
the  ratio  of  free  demand  to  total  dema 
the  'new*  fuel  is  introduced  at  a  rapi 
the  consumption  period  of  a  fuel  is  ex 
free  to  total  demand  will  eventually  s 
which  is  generally  much  smaller  than  1 
TRANSPORTATION  SECTOR* 

As  expected,  the  free  demand  for 
transportation  sector  has  been  very  el 
increase  in  absolute  value  in  the  sixt 
share  of  free  demand  declines.  A  very 
for  the  elasticity  of  total  demand  for 
magnitude  of  elasticity  is  reduced  som 
declining  ratio  of  free  demand  to  tota 

Compared  to  coal,  the  free  demand 
thus,  the  total  demand  for  crude  oil  h 
inelastic  in  this  sector. 

ELECTRICITY  GENERATING  SECTOR1 2 

The  interfuel  substitution  behavi 
generating  sector  is  quite  different  f 


1  see  Figures  IV. 12  to  IV. 16. 

2  see  Tables  IV. 2  to  IV. 4. 
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This  difference  in  behavior  is  explained  primarily  by  the 
economics  of  hydro  generation.  Hydro  generating  projects 
typically  have  large  capital  costs  and  low  operating  costs 
associated  with  them.  Once  the  construction  of  a  hydro 
project  is  complected,  the  economics  of  the  project  dictate 
that  the  plant  be  brought  on-line  at  its  maximum  capacity  as 
scon  as  possible.  This  may  even  result  in  reducing  the 
output  from  fossil  plants  to  below  their  capacity.  It  is  a 
fact  that  the  introduction  of  hydro  plants  in  varying  sizes 
and  at  seemingly  random  time  periods  has  resulted  in  a 
fluctuating  demand  for  fossil  generated  electrical  power.  In 
turn,  this  sector’s  demands  for  crude  oil,  natural  gas,  and 
coal  have  had  corresponding  fluctuations. 

In  analyzing  the  elasticities  of  demand  for  fossil 
fuels  in  this  sector,  one  common  problem  must  be  clarified. 
The  rapid  substitution  of  hydro  power  for  fossil  power  may 
result  in  all  of  the  replacement  demand  of  fossil  plants 
being  captured  by  the  hydro  planus.  In  addition  to  this  loss 
of  fossil  generating  capacity,  the  hydro  generated  power  may 
also  capture  some  of  the  bass  demand  supplied  by  fossil 
plants.  This  will  force  a  reduction  in  the  utilization  rate 
of  the  fossil  capacity.  In  terms  of  the  model  this  simply 
means  that  it  is  possible  for  the  proportion  of  free  demand 
which  is  committed  to  fossil  generation  to  be  negative  under 
such  conditions.  There  can  be  no  meaningful  economic 
interpretation  of  elasticity  of  free  demand  if  this  quantity 
is  negative.  Although  the  utilization  rate  of  the  base 
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It  is  clear  from  the  data  that  the  price 
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would  probably  involve  adjustments  in  their  level  of  energy 
consumption  rather  than  an  immediate  interfuel  substitution, 
A  similar  response  would  be  displayed  by  consumers  in  the 
event  of  a  physical  disruption  of  supplies  of  one  or  more 
particular  fuel  forms. 

The  commitment  liberation  rates  in  the  model  have  been 
assumed  to  be  independent  of  fuel  price  when  in  fact  a 
degree  of  dependence  may  exist.  Consequently,  it  is  believed 
that  the  model  may  be  producing  under -estimates  of  the 
sector  free  demands.  This  is  particularly  true  in  times  of 
rapid  price  changes.  These  biased  estimates  will  be  passed 
on  and  they  will  yield  low  estimates  of  the  sector  demand 
elasticities. 

In  several  cases  the  relative  fuel  prices  are  not  the 
only  determinants  of  the  distribution  multipliers.  What 
then,  are  the  implications  of  the  non- price  factors  in 
regards  to  the  elasticities  of  the  fuel  demands?  Initially 
it  must  be  recognized  that  the  purpose  of  elasticity  studies 
is  to  determine  to  which  influences  a  variable  is  sensitive 
and  to  quantify  the  magnitude  of  this  sensitivity. 

Therefore,  the  elasticity  studies  of  the  sector  demands  are 
net  complete  until  the  variables  representing  TIME  in  each 
equation  are  disaggregated  into  their  basic  components,  and 
the  corresponding  elasticities  estimated.  Although  this 
final  step  has  not  been  taken  here,  the  importance  of  the 
TIME  variables  (and  probably  the  related  elasticities) 
should  not  be  overlooked. 
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FIGURE  IV. 5  -  ELASTICITY  OF  FREE 
COMMERCIAL  SECTOR,  BASE 


DEMAND  FOR  RESIDENTIAL- 
CASE  1945-70. 
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FIGURE  IV. 6  -  ELASTICITY  OF  FREE  DEMAND  FOR  RESIDENTIAL- 
COMMERCIAL  SECTOR,  BASE  CASE  1945-70. 
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FIGURE  IV. 7  -  RATIO  OF  FREE  DEMAND  TC  TOTAL  DEMAND  FOR 
RESIDENTIAL-COMMERCIAL  SECTOR  f  BASE  CASE  1  945-70. 
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FIGURE  IV. 8  -  ELASTICITY  OF 
COMMERCIAL  SECTOR 


OTAL  DEMAND  FOR  RESIDENTIAL- 
BASE  CASE  1945-70. 
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FIGURE  IV. 9 


ELASTICITY  OF  FREE  DEMAND  FOR  INDUSTRIAL 
SECTOR,  BASE  CASE  1945-70. 
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FIGURE  IV,  Tr'  -  RATIO  OF  FREE  DEMAND  TO  TOTAL  DEMAND  FOR 
INDUSTRIAL  SECTOR,  BASE  CASE  1945-70. 
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FIGURE  IV. 11  -  ELASTICITY  OF  TOTAL  DEMAND  FOR  INDUSTRIAL 

SECTOR,  BASE  CASE  1945-70. 
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FIGURE  IV.  12 


ELASTICITY  OF  FREE  DEMAND  FOR  TRANSPORTATION 
SECTOR,  BASE  CASE  1945-70. 
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FIGURE  IV. 13  -  ELASTICITY  OF  FREE  DEMAND  FOR  TRANSPORTATION 

SECTOR ,  RASE  CASE  1945-70. 
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FIGURE  IV. 14  -  RATIO  OF  FREE  DEMAND  10  TOTAL  DEMAND  FOR 
TRANSPORTATION  SECTOR ,  BASE  CASE  1945-70. 
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FIGURE  IV. 15  -  ELASTICITY  OF  TOTAL  DEMAND  FOR  TRANSPORTATION 

SECTOR,  BASE  CASE  1945-7G. 
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FIGURE  IV. 16  -  ELASTICITY  OF  TOTAL  DEMAND  FOR  TRANSPORTATION 

SECTOR ,  BASE  CASE  19U5-79. 
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TABLE  IV. 2 


ELASTICITIES  OF  FREE  DEMAND  FOR  ELECTRICITY  GENERATING 

SECTOR 


YEAR 

COAL 

GAS 

OIL 

1945 

1946 

1947 

1948 

1949 
19  50 

-0,4109 

-0. 9330 

-0.0279 

1951 

-0.3700 

-0, 9019 

-0.0278 

1952 

-0.3165 

-0,8831 

-0.0278 

1953 

-0. 2515 

-0. 8607 

-0.0277 

1954 

-0. 1741 

-0. S350 

-0.0277 

1955 

-0,0836 

-0. 8062 

-0,0277 

1956 

-0.0408 

-0, 7982 

-0.0278 

19  57 

-0.0369 

-0. 7919 

-0.0279 

19  58 

-0,0366 

-0, 7840 

-0.0281 

1959 

-0.0363 

-0. 7745 

-0.0282 

1960 

-0.0361 

-9.7634 

-0.0283 

1961 

-0.0347 

-3.7623 

-0.0288 

1962 

-€.0332 

-0. 761 6 

—  U • 02  92 

196  3 
1964 

-0.3319 

-0. 7593 

-0.0297 

1965 

-0.0543 

-0,  7502 

-0.0308 

1966 

1967 

1968 

1969 
19  70 

-0.2122 

-0. 8630 

-0.0314 

IABLE  IV, 3 


RATIO  OF  FEE  E  DEMAND  TO  TOTAL  DEMAND  FOE  ELECTRICITY 

GENERATING  SECTOR 


YEAR 

1945 

COAL 

GAS 

CIL 

1946 

1947 

1948 

— 

-* 

_  -  — 

19  49 

0.350 

0.441 

0 

.3  36 

19  50 

1951 

0,290 

0.  410 

0 

.334 

19  52 

0. 270 

0.3  75 

0 

,30  9 

1953 

0.  490 

0,  596 

o 

.  52  8 

1954 

0,  6  50 

0,768 

0 

.  694 

1955 

0,710 

0 »  855 

0 

.  752 

1956 

0  •  520 

0.8  75 

0 

.  57  8 

1957 

0.460 

0.  561 

0 

.  476 

1958 

0.660 

0.  795 

0 

.  676 

1959 

0.6  50 

0.769 

e 

.  6  5  6 

1  960 

0,61 0 

0 , 7  74 

o 

.  6  0  €' 

1961 

0.  700 

0 .738 

0 

.  692 

196  2 

0 ,87  0 

0.887 

0 

,  85  8 

1963 

0.9  30 

0.923 

0 

.905 

1964 

— 

- 

— 

1965 

0.9  80 

0.  947 

0 

.  933 

196  6 

— 

- 

— 

1967 

1968 

1,000 

0.868 

0 

.919 

196  9 

■ - 

- 

— 

1979 

— 

- 

- - 

TABLE  IV. 4 
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ELASTICITIES  OF  TOTAL  DEMAND  FOE  ELECTRICITY  GENERATING 

SECTOR 


YEAR 

COAL 

GAS 

OIL 

1945 

1946 

1947 
19  48 

1949 

1950 

-0. 1400 

-0. 4100 

-0.G108 

1951 

-y,  1100 

-0. 3690 

-  .  ,  9  3 

1952 

-0,0900 

-C  .3310 

-0.0086 

1953 

-0.1200 

-C. 5130 

-0.C146 

1954 

-0.1100 

-0.6420 

-0.01 92 

1955 

-0,0600 

-0.6890 

-0,0208 

1956 

-0. 520 

-0.6990 

-0.0161 

1957 

-0.0200 

-0. 4440 

-  0  •  0 1 3  3 

1958 

-o.  :  o 

-0,6230 

-0.01  90 

1959 

-0.0200 

-0. 5950 

-0.0185 

1960 

-0.0200 

-0, 5900 

-C.01 70 

1961 

-C.C2CC 

-v. 5630 

-0.0199 

1962 

-0. 0 300 

-0,6750 

-0.0251 

1963 

1  964 

-0.0  300 

-0. 7010 

-0.0269 

1965 

-0.0500 

-C  .  7 1 1 0 

-0.0288 

1966 

1967 

1968 

1969 

1970 

-0.2100 

-0.7490 

-0,0288 
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C,  PROJECTIONS 


One  majo 
future  projec 
order  to  laeas 
{controllab le 
be  made  on  th 
uncontrollabl 
critical  area 

An  infin 
performed  .by 
under  an  equa 
projections  w 
fuel  prices, 
rates,  fossil 
produced  from 
scenarios. 


r  application  of  quantitative  models  is  making 
tions.  These  projections  may  be  conducted  in 
ure  the  impact  of  changes  in  policy 
)  variables.  On  the  other  hand,  projections  may 
e  basis  of  expected  future  behavior  of 
e  variables  with  the  objective  of  determining 
s  upon  which  government  policy  must  focus. 

ite  number  of  projections  exist  which  could  be 
the  ’Canadian  Interfuel  Substitution  Model’ 

1  number  of  hypotheses.  Some  of  the  interesting 
ould  include  hypothesized  future  scenarios  of 
sector  growth  rates,  commitment  liberation 
capacity  fraction,  the  quantity  of  electricity 
hydro,  and  numerous  combinations  of  these 


However,  since  the  prime  purpose 
demonstrate  possible  uses  of  the  model 
will  be  dealt  with.  These  three  cases 
measure  the  implications  of  two  differ 
scenarios  and  two  different  commitment 
scenarios.  The  justification  for  choos 
includes : 


here  is  simply  to 
,  only  three  key  cases 
will  attempt  to 
ent  future  price 
liberation  rate 
ing  these  scenarios 


-  the  sector  growth  rates  are  not  fundamental  to 
determining  the  extent  or  nature  of  the  interfuel 
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Radical  changes  would  be  required  to  alter  these 
programs  significantly. 

-  Canada  is  currently  in  a  situation  of  high 
uncertainty  with  regard  to  future  expected  fuel  prices 
and  hence  commitment  liberation  rates.  Therefore,  a 
keen  interest  exists  concerning  the  future  reactions  to 
possible  price  movements. 

The  time  span  employed  for  the  projections  is  fifteen 
years.  Since  the  ‘Canadian  Interfuel  Substitution  Model’ 
claims  to  produce  quantitative,  as  well  as  qualitative 
accuracy,  it  was  felt  that  projections  to  the  year  1985 
would  properly  satisfy  these  claims, 

CASE  1 


The  CASE  1  Projection  is  based  on  the  following  future 
scenarios  of  the  exogenous  variables. 

1,  Primary  demand  sector  growth  rates  -  maintained 
equal  to  those  that  existed  for  the  period  1945-70, 

2,  Fossil  capacity  fraction  (FCF)  -  based  on  estimates 
by  the  Department  of  Energy,  Mines  and  Resources1  ’ 2 . 

3,  Quantity  of  electricity  produced  from  hydro  (ZFH)  - 
based  on  estimates  by  the  Department  of  Energy,  Mines 
and  Resources1'3. 

1  Department  of  Energy,  Mines  and  Resources  An  Energy  Policy 
fcr_Canada  Ottawa  Vol.  II,  pg.  292. 

2  see  Figure  IV. 17. 

3  see  Figure  IV. 18. 
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4.  .Fuel  prices  -  assume  that  the  price  of  oil  and  gas 
in  1985  will  triple  the  average  price  that  existed  in 
1 970i . 

Coal  prices  will  rise  significantly  but  net  at  the  same 
rate  as  oil  or  gas  prices.  Electricity  prices  will 
increase  moderately  in  response  to  slightly  higher 
nuclear  generation  costs.  However,  these  rising 
electricity  costs  are  expected  to  reach  a  new  and  very 
stable  plateau  shortly  after  1985, 

5.  Primary  sector  consumer  commitment  liberation  rates 
-  maintained  at  the  1970  values  of  the  Base  Case, 

The  results  of  the  Case  1  projections  are  presented  in 
Figures  IV, 21  to  IV. 27, 

The  fuel  price  scenarios  hypothesized  represent  a  large 
decline  in  the  relative  price  of  electricity  while  the 
relative  prices  of  the  fossil  fuels  remain  effectively 
constant . 

Although  the  responses  of  the  residential-commercial 
sector1 2  and  the  electricity  generating  sector  correspond 
with  the  price  ratio  changes  that  occured,  it  appears  that 
the  industrial  sector3  has  been  dominated  by  the  behavior  of 
the  variable  TIME  (i.e,  those  factors  other  than  relative 


1  see  Figures  IV, 19  and  IV. 2  '  for  the  projected  prices  and 
price  ratios, 

2  see  Figure  IV, 22. 

3  see  Figure  IV, 23, 


97 


fuel 

pri 

ces  that  a if 

ect  id  a 

rk 

et 

sha 

re  b 

eha  vior ) 

in 

the 

distr 

ibu 

tion  multipl 

iers. 

De 

spi 

t  € 

it  s 

lo 

wer  re 

ia  ti 

ve.  p 

rice , 

the  e 

lec 

tricity  mark 

et  sha 

re 

of 

th 

e  in 

du 

st  rial 

sec 

tor 

is 

rapidly 

replaced  by 

natura 

1 

gas 

an 

d  cr 

ud 

e  oil 

in  t 

heir 

res  pe 

cti 

ve  orders  of 

signi 

f  ican 

c  e . 

The 

r 

esult s 

of 

the 

Case  2 

Pro  je 

ctions  further 

discuss 

the 

pr 

oble 

m 

of  the 

non 

~  pri 

ce 

facto 

rs 

in  the  fuel 

select 

io 

n  p 

roc 

ess , 

T  he 

transporta t 

ion  se 

ct 

or 1 

/  a 

nd  p 

ar 

t icula 

rily 

the 

deman 

d  f 

or  crude  oil 

,  has 

be 

en , 

in 

t  he 

P 

ast ,  v 

ery 

insensitive  to  price  changes.  Currently,  and  in  the 
foreseeable  future,  there  are  no  effective  substitutes  for 
crude  oil  in  this  sector. 

The  behavior  of  the  electricity  generating  sector  is 
dependent  on  the  primary  demand  sector  requirements  for 
electricity.  However,  regardless  of  the  behavior  of  the 
primary  demand  sectors,  there  is  a  trend  in  the  electricity 
generating  sector  towards  a  larger  market  share  for  nuclear 
produced  electricity2.  This  trend  exists  by  specification  of 
the  exogenous  model  inputs. 


1  see  Figure  IV. 2 4 

2  see  Figure  IV, 26 
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FIGURE  IV. 18  - 
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FIGURE  IV.  19  -  FUEL  PRICES,  CASE  I  1^45-85. 


PRICE  RATIOS 


101 


-  RELATIVE  PRICES  OF  THE  FUEL  FORMS ,  CASE  I 
1 945-85. 


FIGURE  IV. 20 
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FIGURE  IV. 21  -  ENERGY  DEMAND  BY  THE  PRIMARY  CONSUMING 

SECTORS,  CASE  I  1945-65. 


COM  DEMAND  CQ) 


103 


I 

m 

u 

(Y 


TIME 


(YEARS) 


FIGURE  IV. 22  - 


P  F  5  I D E  N  TI A  L  -  C  0 .1  ME  PC  I A  L 
85. 


DEMAND,  CASE  I  1945- 


INDUSTRIAL  DEMAND  '  (Q 


2  -  0000- 


t 


104 


A  -  Electricity 
V  -  Crude  Oil 
+  -  Coal 
x  -  Natural  Gas 


i » 5000-  - 


x 


X 


1*0000-- 


0 * 5000 


V 


■v 


V 


*7 


0 


ooooixxx^--- 


/X 


X 


X 


X 


X 


X 


XX* ++++++ 


+ 


4- 


4~ 


+4- 


+4-+ 


++ 


in 

O 

in 

O 

in 

O 

in 

O 

m 

NT 

in 

in 

CD 

CO 

N 

N 

CO 

CD 

CD 

CD 

CD 

CD 

01 

Cl 

CD 

CD 

CD 

H 

v4 

H 

H 

H 

H 

•H 

H 

H 

TIME  (YEARS) 


F 


IG11RE 


TV.  2  3 


INDUSTRIAL  DEMAND,  CASE  I  1943-85 


TRANSPORTATION  DEMAND  (Q 


105 


-  TRANSPORTATION  DEMAND,  CASE  I  1945-85. 


FIGURE  IV. 24 
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FIGURE  IV. 25  -  ELECTRICITY  GENERATED  FROM  PRIMARY  FUEL 

FORMS,  CASE  I  19U5-85. 
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CASE  2 

This  case  study  adopts  the  same  scenarios  as  CASE  1 
with  the  exception  of  the  fuel  price  scenarios. 

1.  Fuel  prices  -  this  case  extrapolates  the  19^5—7'' 
fuel  prices  and  completely  ignores  the  unstable  fuel 
price  situation  of  1S73-741. 
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The  results  of  this  projection  are  contained  in  Figures 
IV. 30  to  IV. 38. 
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see  Figures  IV. 28  and  IV, 29. 
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relative  to  ether  factors  as  time  increases.  In  spite  of  the 
fact  that  these  non-price  factors  may  be  powerful,  it  is  not 
logical  to  expect  the  demand  for  a  fuel  to  become 
independent  of  its  price,  particularity  if  fuel  prices 
increase  rapidly  relative  to  other  commodity  prices. 
Unfortunately,  however,  there  is  no  real  basis  for  adjusting 
the  coefficient  of  the  variable  TIME  until  this  variable  has 
been  disaggregated  and  understood  in  its  most  basic  form. 
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FUEL  PRICES,  CASE  II  19U5-85. 
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FIGURE  IV. 30 


-  ENERGY  DEMAND  BY  THE  PRIMARY 
SECTORS,  CASE  II  19U5-65. 
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FIGUFF,  IV.  3  2 


INDUSTRIAL  DEMAND,  CASE  II  1945-85 
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FIGURE  IV. 33 
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FIGURE  IV. 36  - 


:lectricity  generated  from  hydro,  CASE  II 
1945-85. 


MARKE r  SHARES 


TIME 


(YEARS) 


-  PRIMARY  FUEL  MARKET  SHARES  IN  THE 
GENERATING  SECTOR,  CASE  II  1945-85. 
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IV. 38  -  PRIMARY  FUEL  MARKET  SHAKES  FOR  TOTAL 
ENERGY  CONSUMPTION,  CASE  II  1945-85. 
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CASE  3 

The  effects  of  time-varying  commitment  liberation  rates 
fcr  consumers  first  became  a  point  of  interest  during  the 
parameter  estimation  of  the  equations  for  the  distribution 
multipliers.  It  was  discovered  that  time  varying  commitment 
liberation  rates  improved  the  explanatory  power  of  the 
distribution  multipliers,  particularly  in  the  case  of 
industrial  and  transportation  coal  consumers. 

Hence,  this  case  study  examines  the  effect  of  future 
increases  in  the  commitment  liberation  rates  of  crude  oil 
and  natural  gas  consumers.  Since  natural  gas  and  crude  oil 
were  first  introduced  as  significant  sources  of  energy  in 
Canada  in  the  mid  1940* s,  the  consumers  of  the  'new'  fuels 
have  naturally  possessed  very  low  commitment  liberation 
rates  in  the  past.  However,  the  conversion  equipment  of  the 
initial  crude  oil  and  natural  gas  consumers  is  recently 
requiring  replacement  for  the  first  time.  This  is  manifested 
by  increases  in  the  commitment  liberation  rate  (market  entry 
rate)  of  these  consumers.  In  addition,  we  can  expect 
increases  in  the  commitment  liberation  rates  of  oil  and  gas 
as  a  result  of  increases  in  the  relative  prices  of  these 
fuels . 

Therefore,  the  hypothesized  future  scenarios  of  the 
commitment  liberation  rates  for  oil  and  gas  consumers  are 


given  below 
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1,  Commitment  Liberation  rates1 

-  for  the  residential-commercial  sector  ana  the 
industrial  sector,  the  commitment  liberation  ra 
natural  gas  consumers  is  projected  to  increase 
per  year  in  1975  to  20%  per  year  in  1985  while 
commitment  liberation  rate  for  crude  oil  consum 
projected  to  increase  from  10%  per  year  in  1975 
per  year  in  1985.  The  commitment  liberation  rat 
fossil  fuel  consumers  in  the  electricity  genera 
sector  are  projected  to  be  equal  to  those  for  n 
gas  consumers  in  the  residential-commercial  sec 
other  commitment  liberation  rates  are  equal  to 
effect  in  1970  in  the  Base  Case. 

The  other  exogenous  inputs  assume  the  values  th 
assigned  to  them  in  Case  1, 

The  results  of  this  projection  are  presented  in 
IV. 45  to  IV. 53. 

The  consequence  of  increasing  the  consumer  coram 
liberation  rates  in  any  sector  is  to  make  that  secto 
consumption  patterns  more  responsive  to  fuel  price  c 
Obviously,  if  the  commi fluent  liberation  rates  for  al 
consumers  in  a  particular  sector  were  1'0%  per  year, 
the  total  sector  demand  would  respond  instanta neousl 
each  time  period  to  the  fluctuations  in  fuel  prices 
ether  factors.  As  a  result,  the  total  sector  fuel  raa 


1  see  Figures  IV. 32  to  IV. 35 
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shares  would  be  equal  to  the  distribution  multipliers 
without  any  time  delay. 

One  historic  example  of  the  tremendous  impact  that  a 
variable  com mitment  liberation  rate  (regardless  of  what 
causes  the  variation)  can  produce  is  the  sudden  decline  in 
coal  consumption  by  the  transportation  sector  after  1950. 

Although  the  responses  of  the  demand  sectors  in  the 
projections  of  CASE  3  fellow  the  expected  trends,  the  full 
impact  of  the  variable  commitment  liberation  rates  is  once 
again  not  realized  due  to  the  strong  influence  of  the 
variable  TIME  in  the  distribution  multipliers. 

It  is  interesting  to  note,  however,  that  the  results  of 
this  projection  are  similar,  generally,  to  the  results  of 
forecasts  made  by  the  Department  of  Energy,  Mines  and 
Resources1  and  the  Atomic  Energy  of  Canada  Limited  (AECL)2, 
The  projections  of  primary  and  secondary  energy  consumption 
by  fuel  and  secondary  energy  consumption  by  sector  produced 
by  the  'Canadian  Interfuel  Substitution  Model'  are  very 
close  to  the  corresponding  forecasts  made  by  the  Department 
of  Energy,  Mines  and  Resources.  Unfortunately,  further 
comparisons  are  made  difficult  because  of  the  use  of 
different  KWH-BTU  conversion  constants.  Atomic  Energy  of 


1  Department  of  Energy,  Mines  and  Resources  An  Energy  Policy 
foi_Canaua  Ottawa,  1973  pgs,  73-75. 

2  E.  C.  is.  Perryman  "Canadian  Power  Reactor  Pregram  - 
Present  and  Future",  Report  #AECL-4265,  Chalk  River  (1972) 
pg,  1. 
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Canada  Limited  forecasts  that  by  1985  Canada’s  nuclear  power 
capacity  will  be  approximately  14, COO  MWe..  In  CASE  3,  the 
model  projects  that  the  nuclear  power  capacity  of  Canada 
will  be  around  13,000  MWe.  in  1985. 
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FIGURE  IV. 40  -  CONSUMER  COMMI' 
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FIGURE  IV.  41  - 


CONSUMER 


COMMITMENT  LIBERATION  RATES  FOR  THE 


TRANSPORTATION  DEMAND  SECTOR ,  CASE  III  1945-85. 
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FIGURE  IV. 42  -  PRODUCER  COMMITMENT  LIBERATION  RATES  FOR  THE 
ELECTRICITY  GENERATING  SECTOR,  CASE  III  1945-85. 
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FIGURE  IV. 43  -  FUEL  PRICES,  CASE  III  1945-35. 
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FIGURE  IV. 44  - 


RELATIVE  PRICES  OF  THE  FUEL  FORMS,  CASE  III 
1 945-85 . 
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FIGURE  IV. 45  -  ENERGY  DEMAND  BY  THE  PRIMARY  CONSUMING 

SECTORS,  CASE  III  1945-85. 
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FIGURE  TV. 46  -  RESIDENTIAL-COMMERCIAL  DEMAND,  CASE  III  1945- 
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IV. 47  -  INDUSTRIAL  DEWAND,  CASE  III  1945-85 
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FIGIJPE  IV  .  4  9 


ELECTRICITY  GENERATED  FROM  PE 
FORMS,  CASE  III  1945-85. 
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FIGURE  IV. 51  -  PERCENTAGE  OF  TOTAL  ELECTRICAL  GENERATING 
CAPACITY  WHICH  IS  DEPENDENT  ON  NON-NUCLEAR  FUEL,  CASE  III 

1 945-85 . 
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FIGURE  IV. 51  - 


ELECTRICITY  GENERATED  FROM  HYDRO,  CASE  III 
19L5-85. 
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B.  POL  ICY_  1^1  PLICATIONS 

The  broad  interrelated  nature  of  the  energy  industries 
and  their  crucial  importance  in  an  industrial  economy  compel 
the  government  of  Canada  to  make  energy  related  policy  a  top 
priority.  The  ultimate  objective  of  government  energy  policy 
must  be  to  ensure  that  the  energy  requirements  of  Canadians 
are  satisfied  on  a  smooth,  continuous  basis  at  a  minimum 
social  and  economic  cost, 

for  our  purpose  here  it  is  desirable  to  consider  energy 


policy  in  two 

facets  -  demand  policy 

and 

supp 

ly  policy. 

Despite  their 

intimate  relationship. 

t  he 

disc 

ussion  will 

concentrate  on 

energy  demand  policy. 

That 

is , 

what  are  the 

steps  that  may  be  taken  by  the  Federal  Government  which 
directly  affect  energy  consumers  and  satisfy  the  objectives 
of  the  overall  energy  policy? 

1,)  In  1970  more  than  65%  of  Canada’s  total  demand  for 
primary  energy  forms  was  satisfied  by  non -renewable 
resources.  This  dependence  on  non -renewable  energy 
sources  will  inevitably  expand  until  the  time  when 
breeder  reactors  are  technically,  economically  and 
environmentally  feasible.  As  a  consequence  of  the 
depletion  of  these  non- renewable  energy  sources,  there 
will  be  major  interfuel  transitions  in  the  future.  The 
Canadian  government  may  assist  in  making  these 
interfuel  transitions  less  painful  by  making  the 
following  information  readily  available 
a  continuing  casis: 


tc  consumers  on 
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-  What  are  the  short-term  and  long-term  expected  energy 
and  conversion  equipment  costs? 

-  How  secure  are  the  suppliers  of  various  energy  forms? 

-  What  are  the  physical  depreciation  rates  of  various 
conversion  equipment?  In  other  words#  how  long  is  a 
consumer  committed  to  a  particular  fuel  without  an 
unexpected  loss  once  he  makes  necessary  capital 
investment  in  conversion  equipment? 

-  What  are  the  efficiency  specifications  and  relative 
efficiencies  of  conversion  equipment? 

-  In  what  instances  and  at  what  cost  are  capital 
substitutes  for  energy  (i.e.  insulation)  available  to 
consumers? 

The  existence  of  answers  to  such  questions  in  an 
accessible  form  will  aid  consumers  in  formulating 
rational  long-range  stategies  for  energy  utilization. 
This  action  is  seen  as  ’preventive  medicine’  and  would 
be  part icularily  influential  when  aimed  at  consumers 
with  free  demand. 


2.)  The  Economic  Council  of  Canada,  in  its  latest 
Annual  Eeview  (7),  has  recommended  that  the  government 
allow  Canadian  energy  prices  to  rise  to  the  level  of 
their  international  counterparts.  Although  higher 
natural  gas  and  crude  oil  prices  will  achieve  the 
objectives  of  conservation  and#  where  possible,  force 
the  interfuel  substitution  process#  such  an 
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indiscriminate  adjustment  in  prices  may  net  be 
consistent  with  the  goals  of  satisfying  energy  demand 
at  minimum  social  and  economic  cost.  That  is,  the  same 
goals  and  objectives  may  be  better  achieved  through  the 
identification  of  four  basic  energy  consuming  groups. 
The  base  demand  and  the  free  demand  energy  consumers 
constitute  two  distinct  groups  of  consumers  while  the 
efficient  and  inefficient  energy  consumers  form  two 
more  distinct  groups  of  consumers.  It  is  suggested  that 
significant  changes  in  Canada’s  total  energy 
consumption  patterns  can  arise  from  policies  which 
single  out  and  influence  the  free  demand  consumers  and 
the  inefficient  consumers,  A  government  policy  which 
produces  large  energy  cost  increases  for  efficient  base 
demand  consumers  is  not  justified.  Such  a  policy  will 
not  alter  the  consumption  patterns  of  these  consumers 
without  drastic  side  effects.  On  the  other  hand, 
inefficient  and  free  demand  consumers  will  exhibit  very 
elastic  behavior  with  respect  to  energy  costs  and 
alterations  in  consumption  patterns  will  be  made  with 
relative  ease. 


CHAETEE  V 


CONCLUSIONS  AND  FUETHES  RESEARCH 

In  this  study  a  preliminary  model  of  the  interfuel 
substitution  process  in  Canada  has  been  developed.  The 
model,  as  it  now  exists,  is  referred  to  as  a  preliminary 
ncdel  because  the  developments,  to  date,  simply  identify  the 
major  factors  or  components  in  the  interfuel  substitution 
process  without  really  analyzing  these  components  in  depth. 

These  critical  factors  in  the  inter fuel  substitution 
process  are  described  below  along  with  guidelines  for 
further  research. 

CCN.  MITHEKl  LIBERATION  RAIES  -  The  broad  definition  that  has 
been  employed  for  the  consumer  commitment  liberation  rates 
makes  further  analysis  of  these  rates  very  cc&plex.  In  any 
event,  the  consumer  commitment  liberation  rates  have  been 
clearly  identified  in  this  study  as  the  major  factors  that 
are  capable  of  altering  consumer  flexibility  and 
responsiveness.  It  is  expected  that  in  addition  to  being 
dependent  on  the  relative  prices  of  energy  fuels,  these 
ccaiitmert  liberation  rates  are  also  affected  by 
technological  considerations  (i.e.  inventions  and 
innovations  in  conversion  equipment)  and  tastes.  However, 
further  research  of  the  consumer  commitment  liberation  rates 
should  probably  proceed  on  a  disaggregated,  sector-by-sector 
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, 


. 


.  . 

basis  and  even  include  seme  case  studies. 


DISTRIBUTION  EULTIPIIESS  AND  RELATIVE  EUFL  PRICES  -  During 
the  development  cf  the  distribution  multipliers  it  became 
clear  that  fuel  prices  as  well  as  some  non-price  factors 
entered  Into  the  consumer's  decision  process.  A  significant 
result  of  this  study  was  the  finding  that  relative  fuel 
prices,  rather  that  absolute  fuel  prices,  represent  the  fuel 
price  aspects  in  the  decisions  cf  free  demand  allocation.  In 
terms  cf  a  nationally  aggregated  model,  the  role  of  fuel 
prices  in  the  energy  consumer  decision  process  is  not 
considered  a  major  area  for  further  research. 

DISTRIBUTION  MULTIPLIERS  AND  NON-PRICE  FACTORS  -  Contrary  to 
the  assumptions  made  by  Baughman  (1),  it  has  been  found  in 
this  study  that  the  behavior  of  the  distribution  multipliers 
is  not  adeguately  described  by  the  fuel  prices  alone. 
Although  the  additional  explanatory  variables  in  this  study 
were  represented  simply  by  a  variable  TINE,  it  is  believed 
that  this  variable  is  a  manifestation  cf  other  underlying 
factors.  It  is  difficult  tc  speculate  on  the  fundamental 
determinants  of  these  non- price  factors.  This  is  an 
important  area  for  further  research  and  it  is  expected  that 
the  results  will  be  different  for  each  distribution 
multiplier.  Another  area  of  interest,  which  was  not  covered 
in  this  study,  relates  to  the  time  trends  and  interactions 
cf  the  relative  weightings  cf  the  price  and  ncr-price 
factors  in  the  distribution  multipliers. 


. 


. 


, 


UTILIZATION  BAIE  -  Khazzccm  (11)  has  pointed  cut  that 
consumers  may  employ  variable  utilization  rates,  for 
example,  a  consumer  with  replacement  demand  may  decide, 
after  entering  the  marketplace,  not  no  satisfy  this  demand. 
He  may  prefer  to  leave  his  conversion  equipment  temporarily 
inactive.  Unfortunately,  this  option  does  not  exist  in  the 
present  structure  of  the  "Canadian  Interfuel  Substitution 
Model".  However,  this  factor  could  easily  be  incorporated 
into  the  model  structure  if  the  determinants  cf  variable 
utilization  rates  were  fully  understood. 

E  BICE  ELASTICITY  OE  EISTBIEUTICK  MULTIPLIERS  -  The  discovery 
that  the  price  elasticities  of  the  fuel  demands  vary 
inversely  with  the  fuel  market  shares  is  not  completely 
understood.  To  aid  in  understanding  this  phenomenon  it  would 
be  helpful  to  find  other  examples  where  the  price 
elasticities  cf  commodity  substitutes  are  dependent  on  their 
relative  market  shares. 


Substantial  credit  must  be  given  to  the  methodology  of 
large-scale  mathematical  modeling  for  assisting  in  the 
qualitative  and  quantitative  identification  of  the  major 
factors  in  the  interfuel  substitution  process  and  their 
interrelationships.  In  addition,  the  model  facilitates 
continuing  work  within  a  uniform  framework.  Seme  natural 
extensions  and  recommendations  for  future  investigations  of 
the  "Canadian  Interfuel  Substitution  Model"  are: 


. 


1 

(i)  Development  cf  a  price  determining  dynamic 

ID  C  d  €  1  • 

{ii)  Development  cf  a  dynamic  supply  model. 

Ecth  of  {i)  and  (ii)  above  could  be  easily  adopted  from 
work  cf  Baughman  { 1 ) . 

(iii)  Further  research  into  the  determinants  of 
total  and  per  capita  sector  demands. 

(iv)  Disaggregation  of  the  ’’Canadian  Interfuel 
Substitution  Model”  to  include  consideration  of 
regional  phenomenon  in  Canada. 

(v)  An  important  issue  which  has  not  been  properly 
addressed  in  the  past  relates  to  the  matter  of  BTU 
conversion  constants.  Since  these  conversion 
factors  are  dependent  upon  conversion  equipment 
efficiency  rates,  they  should  be  regularly 
reviewed  and  updated.  This  seems  particularly  true 
for  the  KWE-BTU  conversion  constants  which  now 
assume  two  distinct  values,  depending  on  the  time 
period  in  question. 


. 


. 

. 

. 
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This  appendix  contains  a  presentation  and  discussion  of 
the  data  upon  which  the  'Canadian  Interfuel  Substitution 
Model'  is  based.  Although  the  data  were  converted  to  the 
ccumon  eguivilents  of  British  Thermal  Units  (BTU's)  in  the 
model,  only  the  data  in  its  original  units  will  be  presented 
here.  However,  the  conversion  constants  used  to  convert  the 
data  from  its  natural  units  to  BTU's  will  be  given. 

The  sequence  in  which  the  data  base  will  be  presented 
and  discussed  will  proceed  frcm  the  fuel  price  data  to  the 
sector-b y-sectcr  demand  data  for  each  fuel  forn  and  finally 
tc  the  conversion  constants. 


T  U  E I  PBICE  DATA 

'The  'Canadian  Interfuel  Substitution  Model'  was 
constructed  at  a  high  level  cf  aggregation.  It  was  felt  that 
the  fuel  prices  most  easily  represented  and  corresponding  to 
this  level  of  aggregation  would  be  the  fuel  prices  that 
exist  at  the  site  cf  their  respective  production.  That  is, 
the  prices  of  oil  and  gas  would  correspond  to  the  'well¬ 
head'  prices  and  the  price  cf  coal  would  be  the  price  paid 
for  it  at  the  mine  site.  Unfortunately,  the  data  were  not 
available  to  calculate  the  price  of  electricity  in  this  same 
manner.  As  a  result  the  price  of  electricity  actually 
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. 


. 


represents  the  average  price  paid  by  consumers  in  Canada. 


Tables  A .  1  to  A.  5  contain  the  data  required  to 
calculate  the  particular  fuel  prices.  In  order  to  translate 
the  fuel  prices  into  constant  dollars  the  wholesale  price 
index  fcr  non-metalic  minerals  (1935-39=100)  was  employed. 
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TABLE  A » 1 
C  C  A I  PRICES 


TOTAL  CANADIAN 


EECEOCIION  OF 

VAIOE  OF 

price 

PRICE 

CCAL 

TOTAL  PRODUCTION 

PER  ION 

PER  TON 

YEAR 

(S.10NS  x  106) 

(DOLLARS  x  10^} 

( CDF  PENT  $) 

(1935-39  $) 

1  9  *4  5 

16.507 

67.588 

4.0  94 

3.517 

1946 

17.812 

75.82 

4.256 

3.625 

1947 

15. 869 

77.475 

4.8  82 

3.781 

19  48 

18.45 

106.684 

5.782 

3.  834 

1949 

19.12 

110.915 

5 . 8  C  0 

3.  6  64 

1950 

19.139 

110.14 

5.754 

3.494 

1951 

18.587 

109.039 

5.  866 

3.454 

19  52 

17.579 

111.026 

6.315 

3.631 

1953 

15.901 

102.722 

6.460 

3.  6  51 

1954 

14.914 

96.6 

6.477 

3.659 

1955 

14.819 

93.579 

6.314 

3.604 

19  56 

14.916 

95.35 

6.392 

3.53  5 

1957 

12.96 

90.252 

6.963 

3.  678 

1958 

11.627 

8C  .059 

6.885 

3.6  52 

1959 

10.514 

74.089 

7.046 

3.778 

1960 

10.776 

74.879 

6.948 

3.  743 

1961 

10.336 

70.181 

6. 789 

3.66  6 

1962 

10.217 

6  9.2 

6.773 

3.5  81 

1  963 

10.452 

72.052 

6.893 

3.  637 

1964 

11.219 

73.013 

6.5C7 

3.  409 

1965 

11.5 

76.295 

6.634 

3.462 

1966 

11.18 

81.801 

7.316 

3.777 

1967 

11.141 

56.5 

5.071 

2.545 

1968 

1G.989 

53.97 

4.911 

2.384 

1969 

10.672 

50.578 

4.739 

2.256 

1970 

16 .604 

86.067 

5.183 

2.401 

1971 

18.432 

1 2  1.727 

6.6  04 

2.925 

1972 

20.949 

154. 151 

7.356 

3.  150 

1973 

21.960 

176.979 

8.059 

3.  156 

1AELE  A  «  2 
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NATURAL  GAS  PRICES 


TOTAL  CANADIAN 


PRODUCTION  OF 

V  A  ID  E 

OF 

PRICE  PER 

PRICE 

PER 

NATURAL 

GAS 

TOTAL  PRODUCTION 

1,000  ft. 3 

1  ,000 

ft.3 

¥E  AR  (CD. FT. 

x  109) 

{DOLLARS 

x  10«*) 

(CURRENT  $) 

(1935- 

3  9  $) 

1945 

43.412 

12.31 

0.254 

0.218 

1946 

47  .9 

12.165 

0.253 

0.216 

1947 

52.657 

13.43 

0.255 

0.  197 

1948 

58.603 

15.633 

0.266 

0.  1  76 

19  49 

60.457 

11.62 

0.192 

0.  121 

195G 

67.822 

6.433 

0.094 

0.057 

1951 

79.461 

7.159 

0.090 

0.053 

1952 

88.686 

9.518 

0.  107 

0.061 

1953 

ICO. 966 

10.877 

0.  1  07 

0.060 

1954 

120.735 

12.482 

0. 1  C  3 

0.058 

1955 

150.772 

15.099 

0.  ICO 

0.057 

1956 

169.153 

16.85 

0.099 

0.  055 

1957 

220  .007 

20.863 

0.095 

0.050 

1958 

337.804 

32.058 

0  .  G94 

0.  050 

1959 

417.335 

39.609 

0.094 

0.  050 

I960 

522.972 

52.197 

0.099 

0.053 

1961 

634.131 

68.422 

0.107 

0.058 

1962 

894.672 

108.641 

0.131 

0.064 

196  3 

993  .  38  8 

150.469 

0.151 

0.  079 

1964 

1134.211 

145.658 

0.  128 

0.067 

1965 

1236.798 

158.938 

0.  128 

0.067 

1966 

1341.831 

178.1 84 

0.132 

0.068 

1967 

147  1 .735 

197.983 

0.134 

0.067 

1966 

1696.68 

225.648 

0.133 

0.064 

1969 

1977.838 

262.856 

0.132 

0.063 

197  0 

2276.579 

315.1 

0.138 

0.064 

1971 

2499.024 

342.549 

0.137 

0.061 

1972 

2913.824 

388.905 

0.  133 

0.057 

1973 

5152.410 

482.155 

0.153 

0.060 

TABLE  A. 3 


CRUDE  CIL  PRICES 


TOTAL  CANADIAN 


PRODUCTION  OF 

VALUE  CF 

PRICE 

PRICE 

CRUDE  CIL 

TOTAL  PRODUCTION 

PER  BARREL 

PEP  BARREL 

7  E  AB  (ELLS,  x  10*) 

(LOILAF.S  x  106) 

(CUFF.  ENT  $) 

(1  935-39  $ 

1  9  U5 

9.483 

13.632 

1 ,6  06 

1 . 380 

1946 

7.586 

14.989 

1.975 

1.683 

19  47 

7.692 

19.576 

2.544 

1.971 

19  4  8 

1  1  .896 

37.,  419 

3.145 

2,  085 

19  49 

21.011 

61.118 

2.5C8 

1 . 837 

1950 

28.646 

84 . 6 

2.953 

1 . 793 

1951 

47.535 

116.7 

2.455 

1.445 

1952 

61.161 

143. 

2.3  38 

1.  344 

1953 

80.78 

200.6 

2.4  83 

1 . 403 

1954 

95.96 

243.9 

2,  54  1 

1 . 435 

1955 

129.26 

3C5.6 

2.364 

1  .  349 

1956 

171  .758 

4  06.6 

2.367 

1 . 309 

1957 

181.127 

4  53.6 

2. 5G4 

1 . 322 

1958 

164.679 

39  8.7 

2.421 

1 .2  84 

1959 

183.601 

422.1 

2.299 

1.232 

I9  60 

189.515 

42  2,9 

2.231 

1 . 202 

1961 

22C.816 

487.6 

2.2C8 

1 .192 

1962 

244.195 

5  52.4 

2 

1.196 

1963 

258. 162 

€15.2 

2.382 

1 . 257 

1964 

275.417 

6  7  6.3 

2.  455 

1 .  286 

196  5 

2  9  2  .  j  2  2 

719. 

2.459 

1 . 283 

1966 

320.549 

789.3 

2,462 

1.271 

1967 

251.288 

866  . 

2.465 

1.237 

19  6  8 

373.492 

9  37.4 

2.470 

1 .199 

19  69 

41C.99 

1014.6 

2.468 

1.175 

1970 

461.18 

1156.5 

2,5  C7 

1.162 

1971 

492.739 

1356.9 

2.754 

1.220 

1972 

562.31 

1570.8 

2.793 

1.196 

1973 

649.868 

2246, 149 

3,456 

1.353 
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TABLE  A. 4 

ELECTRICITY  PRICES 


TOTAL  CONSUMPTION 

TOTAL  REVENUE 

PRICE 

PRICE 

CF  ELECTRICITY 

FROM  CUSTOMERS 

PER 

PER 

IN  CAN ALA 

IN  CANADA 

KWH. 

KWII. 

YEAR  (KWH,  .x  109) 

(DOLLARS  x  109) 

(CURRENT  $) 

(1935-39  $) 

15  45 

34.7 

0,215 

0 .006 

0.005 

1946 

37,5 

0.266 

0.0  C7 

0.006 

1547 

41. 

0.244 

0.0  05 

0.004 

1548 

45.1 

0,2  5  7 

0.CC5 

0.003 

1949 

48.9 

0.28 

0.CC5 

0.003 

195  0 

53.459 

0.324 

0.0  06 

0.003 

1951 

59.409 

0.37  5 

0.0C6 

0.003 

1952 

6  4  .  C  2  3 

0.415 

0 .0  C6 

0.003 

195  3 

68.057 

0.469 

0. 0  C6 

0.003 

1954 

71.94 

C  .506 

0.0C7 

0.003 

155  5 

77  .946 

0.549 

0 .007 

0.004 

1956 

83.59 

0.597 

0.0  07 

0.003 

1957 

87.115 

0.64 

0. 0C7 

0.003 

1958 

93.685 

0.692 

0.007 

0.003 

1659 

ICC. 59 

0.756 

0,007 

0.004 

1960 

109.302 

0.805 

0.007 

0.003 

1961 

110.95 

0.85  9 

0 .007 

0.004 

19  62 

1 16. 135 

C.9C8 

0.0  07 

0.004 

1963 

121  .51 

0.966 

0.007 

0.004 

1964 

133.649 

0.963 

0.007 

0.003 

196  5 

144.165 

1.043 

0.0  07 

0.003 

196  6 

156.956 

1  .131 

0.007 

0.003 

16  67 

165.812 

1,248 

0 . 0  07 

0.003 

1968 

176,841 

1  .386 

0.007 

0.003 

1569 

189.522 

1 .53 

0  0  C  8 

0.003 

1 97C 

202 . 337 

1.716 

0.  CC8 

0.003 

15  71 

212.530 

1  .878 

0 . 0  C8 

0.003 

1972 

231.557 

2.065 

0. 008 

0.00  3 

SABLE  A ,5 


WHOLESALE  PRICE  INDEX 


WHOLESALE  PRICE 
INDEX  FOE 
NON-MIL.  MINERALS 


YEAR 

(1935-39  = 

1945 

1.164 

1946 

1  .  174 

1947 

7.291 

15  4  8 

1.508 

1949 

7 .  583 

1  95  C 

1  .647 

1  951 

1  .,698 

1952 

1  .739 

195  3 

1 .7  69 

1  954 

1.77 

1  Q  E  c 

1  y  J  J 

1.752 

1956 

1  .808 

1  957 

1 . 893 

1958 

1.885 

1959 

1  .  865 

I960 

1.656 

1961 

1 . 852 

1962 

1  .891 

1  963 

1 . 895 

1  964 

1  .909 

1  965 

1.916 

1  96  6 

1.937 

1  967 

1.992 

1968 

2.06 

1569 

2.  1 

1  97  C 

2.158 

1971 

2.258 

1972 

2.3  36 

1973 

2.554 
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The  sector-by-sector  demand  statistics  were  difficult 
tc  obtain  cn  a  uniform  basis  because  of  the  changing 
classifications  in  the  data.  As  a  result,  many  adjustments 
were  required  in  the  published  data  in  order  that  it  would 
conform  tc  the  data  requirements  of  the  model.  These 
adjustments  are  outlined  in  the  following  discussions. 

EESIDEHTIAL-CCMHERCI&L  DEMAND  DATA 

The  residential-commercial  demand  data  for  natural  gas 
(1S50-70)  is  equal  tc  the  sum  of  the  residential  demand  data 
and  the  commercial  demand  data  (28) .  Since  the  data  for  the 
1945-50  time  period  were  net  available,  the  statistics  were 
obtained  by  simple  extrapolation  of  the  1950-7C  statistics. 

The  demand  data  for  electricity  were  obtained  in  a 
mar.ner  similar  to  that  for  natural  gas.  The  residential- 
ccmmercial  demand  data  for  electricity  (1950-7C)  is  the  sum 
cf  the  residential  and  commercial  demand  data  (30).  The  data 
fer  1945-50  were  obtained  by  extrapolating  the  later  data. 

The  data  for  residential-commercial  coal  demand  were 
available  (27)  in  the  desired  form.  It  was  necessary, 
however,  tc  extrapolate  the  1949-70  data  in  order  to  obtain 
the  statistics  for  the  period  1545-49. 

Complex  problems  were  encountered  with  the  data  for  the 
residential-commercial  demand  fer  crude  oil.  The  data  for 
the  period  1945-62  were  available  (29)  in  the  desired  form. 
The  data  for  1565-67  were  obtained  by  adding  the  separate 
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classif icat lens  of  residential  demand  and  commercial  demand 
f c  i  coal.  Then  the  values  for  the  period  1  963-64  were 
estimated  by  interpolation  and  the  values  for  the  period 
1966-70  were  estimated  by  extrapolation* 

Table  A, 6  contains  the  data  as  it  was  used  in  the 
model. 

IMUSTBIAL  CEKANC  EATA 

There  were  greater  problems  encountered  in  obtaining 
the  data  for  the  industrial  demand  sector  in  its  desired 
fcim  than  any  ether  demand  sector.  This  was  particularly 
true  for  the  coal  and  crude  oil  consumption  data. 

The  industrial  demand  data  for  natural  gas  (28)  and 
electricity  (30)  were  available  in  the  proper  form  in  the 
original  sources.  It  was  necessary  in  both  cases,  however, 
tc  extrapolate  the  1950-70  data  in  order  to  obtain 
statistics  for  the  period  1945-50. 

The  industrial  demand  data  (29)  for  crude  oil  for  the 
period  1Q65-70  were  acquired  by  subtracting  the  commercial 
demand  fer  crude  oil  from  the  industrial-commercial  demand 
for  crude  oil  (29)  .  The  data  for  the  period  1945-65  was  then 
estimated  extrapolating  the  1965-70  data.  The  extrapolation 
was  based,  tc  a  large  extent,  cn  the  demand  trends  that 
existed  in  other  sectors  for  crude  oil.  It  was  also  based  on 
the  recognition  that  crude  oil  is  a  relatively  ’new'  fuel  in 
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For  the  period  1963-70,  the  industrial  demand  for  coal 
exists  in  the  required  form  (27) .  Before  this  time  period 
electricity  generating  demand  for  coal  and  the  industrial 
demand  for  ccal  were  classified  in  one  category.  Therefore, 
the  electric  utilities  demand  (27)  for  coal  was  subtracted 
frcm  the  •industrial-electricity*  (27)  demand  fcr  coal  to 
obtain  the  required  statistics  fcr  the  period  1949-62.  The 
unavailable  statistics  fcr  the  period  1945-49  were  obtained 
through  extrapolation  of  the  later  data. 


Table  A. 7  contains  the  industrial  demand  data  as  it  was 
used  in  the  medal. 


TRANSPORT ATION  DEMAND  DATA 

The  transportation  demand  for  coal  existed  almost 
solely  within  the  railway  sub-sector.  Consequently,  the 
railway  demand  data  for  coal  (27)  was  collected. 

The  transportation  demand  for  crude  oil  (29)  was 
assumed  tc  consist  cf  the  sum  of  the  road,  rail,  marine,  and 
aviation  demands  for  this  fuel.  Unfortunately,  these  data 
were  available  in  the  required  form  only  for  the  period 
1965-70.  The  remainder  cf  the  statistics  were  obtained  by 
extrapolation  of  the  1965-70  data. 

Since  the  transportation  demand  for  electricity  and 
natural  gas  is  negligible,  nc  attempt  was  made  to  collect 


this  data. 
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Table  A*  8  contains  the  transportation  dens  and  data  as  it 
was  used  in  the  model. 


El  ECTEICI 1Y  GENERATING  DEMAND  DATA 

Electricity  is  usually  considered  to  be  generated  by 
cne  of  three  broad  categories:  public  utilities,  private 
utilities  and  industrial  establishments. 

The  only  available  data  for  the  electricity  generating 
demand  sector,  with  the  exception  of  hydro  generation,  was 
in  the  time  period  1956-70,  In  the  case  of  hydro  generation, 
the  data  were  available  for  the  period  1949-70*  As  a  result, 
estimates  of  the  required  data  were  obtained  by 
extrapolation  wherever  necessary. 

All  the  necessary  data  for  this  sector  are  available 
from  one  source  (30)  and  are  contained  in  Table  A. 9. 

ETC  CONVERSION  CONSTANTS 

The  ETU  conversion  factors  are  functions  which 
translate  the  energy  quantities  from  their  natural  units  to 
Eritish  Thermal  Units  (ETO's) ,  With  the  exception  of 
electricity,  these  functions  have  been  represented  by  all 
researchers  since  1 9 b 5  as  time  invariant  constants. 

However,  it  has  been  recognized  (10)  that  the 
efficiency  cf  the  electrical  generation  process  has  been 
steadily  improving.  As  a  result,  two  conversion  constants 
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have  arisen.  They  are  based  cn  the  estimated  conversion 
efficiencies  cf  fossil  fuel  generating  plants  in  the  years 
1955  and  1968.  The  'Canadian  Interfuel  Substitution  Model' 
utilizes  a  K8H-B1U  conversion  constant  which  lies  between 
the  two  accepted  values  but  is  much  closer  to  the  1968 
estimate . 

The  constants,  as  used  for  the  model,  are  listed  in 
Table  A. 1C  and  the  units  correspond  to  the  natural  units  of 
the  data  as  collected* 
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TABLE  A »  6 

RESIDENTIAL-COMMERCIAL  DEMAND  FOE  E  NEE  GY 


NATURAL 

CRUDE 

CCAL 

GAS 

OIL 

ELECT  RICH 

YEAR  {S.TONS  x  106) 

(CU.FT.x109) 

(BELS. xl 0* ) 

(K¥H„ xIO9) 

1945 

1C.  5 

18. 

4.711 

5.  9 

1946 

10.3  1 

2  2  » 

9.969 

6.9 

1947 

1C.1 

26. 

16.271 

8.  1 

1948 

5.8 

31  . 

17.034 

9.6 

1949 

5.259 

36. 

18.731 

10.9 

1950 

5.06 

40.004 

24.669 

1 2. 685 

1551 

8.55 

43.048 

29.786 

13.939 

1952 

7.144 

43.328 

34.863 

16. 005 

19  53 

6.81 

46.391 

38.586 

17.439 

1954 

6.554 

56.864 

46.808 

19.618 

19  5  5 

6.724 

68.591 

52.862 

2 1. 992 

1956 

6.3  3 

77.937 

61.277 

24.  16 

1957 

5.513 

92.218 

63.17 

26.064 

19  58 

5.225 

112.94 

68.108 

27. 865 

1559 

4.747 

142.683 

74.003 

31.023 

1560 

4.2  15 

161.255 

77.375 

33.029 

1561 

3.997 

175.937 

81.341 

36.391 

1962 

3.617 

198.594 

87.942 

39. 816 

1563 

3.336 

216.219 

101. 

43. 916 

1964 

3.03 

247. 101 

1  14  . 

43.702 

196  5 

2.609 

288.348 

130.032 

53.597 

1966 

2.187 

31 C. 983 

1 32.603 

61.218 

1567 

1.907 

335.775 

137.206 

6  6.  2  4 

1968 

1.533 

359.248 

141.7 

71.932 

1969 

1.306 

399.885 

1  46. 

78. 935 

19  7  C 

1.016 

434.61 

1  50. 

37.499 

TABLE  A. 7 


INDUSTRIAL  DEMAND  FOR  ENERGY 


NATURAL 

CRUDE 

COAL 

GAS 

OIL 

FLECTFIC 

YEAR  (S . TG  N  Sx  1 06) 

(CU.FT.x10g) 

(BBLS. xl 06 ) 

jKNH.xlO 

1945 

1 1.9 

3  . 

19.113 

27. 

1946 

11.65 

6  . 

18.780 

28.5 

1947 

11.45 

10. 

15.805 

30.2 

19  48 

11.3 

13. 

13.679 

32. 

1949 

11.421 

i  c; 

i  —  » 

19. 035 

34. 

1950 

11.601 

17.897 

18.750 

35.774 

1951 

11.714 

21.869 

2  1„ 36  8 

39.691 

1952 

1  1.357 

22.677 

21.311 

42.01 

19  53 

1  1.986 

24.284 

21.239 

44.184 

1954 

10.504 

30.275 

24.502 

45.523 

1955 

11.322 

48.699 

28. 877 

48.634 

1956 

1  1.564 

65.631 

34.907 

51.161 

1957 

10.581 

75.385 

35.307 

52.583 

1958 

9.562 

90.08 

36. 287 

56.876 

1959 

8.782 

1 39.567 

40.464 

59.736 

1960 

8.386 

164.234 

43.459 

66.353 

1961 

7.958 

194.803 

48. 487 

64.354 

19  62 

7.673 

213.468 

52.166 

65.498 

1963 

7.715 

235.379 

58.015 

66.579 

1964 

8.081 

257.402 

59. 523 

72.694 

1965 

8.279 

284.669 

62. 534 

77.28 

1966 

7.543 

324.532 

67  .854 

81 . 546 

1967 

7.381 

359.331 

74.654 

84.539 

1968 

6.56 

406.578 

86. 593 

86.068 

1969 

5.442 

444.829 

86.917 

94.435 

1970 

3.7  18 

399.057 

108.638 

97.02 

TABLE  A. 8 


Tfi AN SPORT AT ION  DEMAND  FOR  ENERGY 


CRUDE 


COAL 

OIL 

YEAR  (S.TONSxlO6) 

(BBLSxlO6) 

1845 

12.084 

38.226 

1846 

11.632 

37.560 

1847 

12.33 

39.513 

1948 

12.422 

45.365 

1849 

12.855 

49.491 

1950 

10. 95 3 

56. 251 

195  1 

10.712 

70.145 

1952 

9.52 

76.801 

1  553 

8.002 

83.564 

195  4 

6.33 

91.579 

195  5 

5.852 

103.957 

1956 

5.324 

127. 992 

1957 

3.109 

134.855 

1958 

1.658 

13  3.  79  5 

1859 

0.874 

146.443 

1560 

0.638 

152. 630 

1561 

0,466 

159.201 

1962 

0.375 

163.951 

1563 

0.307 

169.801 

1964 

0.  271 

162. 450 

1  96  5 

0.22  3 

154.061 

1  966 

0.212 

171.538 

1  567 

0.162 

166.324 

1968 

0.132 

172.963 

1569 

0.132 

187.  123 

1570 

0.138 

197.215 

* 
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TABLE  A. 9 

ELECTRICITY  EEOEUCED  EY  PRIMARY  ENERGY  SOURCES  {KWH.  x  10^) 


YEAR 

CCAL 

NATURAL 

GAS 

CRUDE 

OIL 

NUCLEAR 

HYDRO 

1945 

0.603 

0.026 

0.069 

0. 

34. 

1946 

C.593 

C.031 

C„  C74 

0. 

3  6.8 

1  947 

G  .  6  6  2 

0.043 

0.094 

0. 

40.2 

19  48 

0.729 

C  .  056 

0.113 

0. 

4  4.2 

1949 

4.839 

0.441 

0.839 

0. 

42.779 

1950 

5.245 

0.590 

0.999 

0. 

46.624 

1951 

4.733 

0.6  94 

1. 025 

0. 

51.955 

1952 

4.899 

0.8  92 

1 .207 

0. 

57.024 

1953 

6.102 

1 .  401 

1.627 

0. 

58.926 

1954 

5.697 

1.701 

1.769 

0. 

62.572 

1  955 

4.992 

1.857 

1 .617 

0. 

69.478 

1  956 

1.173 

0.56  8 

0.438 

0. 

81,408 

1957 

2.0  74 

1.081 

0.585 

0. 

82.373 

19  53 

1.468 

1.276 

0.431 

0. 

9C  .50  9 

1959 

1  .272 

1.783 

0.4  94 

0. 

97.04 

1960 

1.403 

1.565 

0.4  49 

0. 

1 05.883 

1961 

2.926 

3.056 

1  „  048 

0. 

102.919 

1  962 

7.251 

2.394 

1.411 

0.027 

104.051 

1963 

12.182 

3.958 

1 .432 

0.  104 

102.832 

1  964 

14.532 

3.875 

2  .  C  3  8 

0.169 

113,334 

1  965 

13.371 

5.572 

3. 016 

0.  141 

117.063 

19  66 

17.821 

6.875 

3.236 

0.188 

129.834 

1  967 

21.538 

6.250 

5.111 

0.164 

132.747 

1968 

2  5 . 7  G4 

7.233 

7.958 

0.971 

134.973 

1969 

27.870 

4.794 

7.048 

0.560 

149.247 

1970 

22.975 

5.768 

5. 8C3 

1.080 

156.709 
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TABLE  A  » 1 0 

ETO  CC  N  V  F  F  S  JC  N  CONSTANTS* 


CCAL 

1  x  1  0  6  short  tens  of  coal  =  .026  Q 


NATURAL  GAS 

1  x  10q  cubic  feet  cf  natural  gas  =  .00107  Q 


CRUDE  CII 

1  x  106  barrels  cf  crude  cil  =  „0058  Q 


IIICTEICITY 

1  x  109  kwh«  cf  electricity  -  .CG94  Q 


*  1  .0  C  =  1  x  10*5  ETU 


. 


, 


. 
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APPENDIX  3 


PASS  BET  EE  ESTIMATION  OF  DISTRIBUTION  MULTIPLIERS 

This  appendix  presents  the  methodology  and  final 
results  cf  the  parameter  estimation  for  the  distribution 
multiplier  equations . 

The  data  which  was  required  for  the  parameter 
estimation  consisted  cf  the  proportion  of  free  demand  which 
was  satisfied  by  each  fuel  in  each  time  period.  In  order  to 
acquire  this  data  cn  a  sector-by-sector  basis,  it  was 
necessary  to  utilize  the  theoretical  structure  cf  the 
•Canadian  Interfuel  Substitution  Model*  in  conjunction  with 
the  assumptions  made  relating  to  consumer  commitment 
liberation  rates.  The  following  equations  present  the 
development  of  an  equation  for  calculating  the  data  for  the 
actual  distribution  multipliers.  This  procedure  is  described 
for  the  residential-commercial  sector  but  it  applies  equally 
to  all  the  demand  sectors. 

The  methodology  begins  with  the  model  equation  for 
residential-commercial  demand  for  fuel  m. 

ECEm.K  =  BCD  m . J  +  BT*[ -ECDmDE. JK 

+  ECEDa.J*  RCMSD.JK]  (B . 1 ) 

In  differential  equation  form  this  equation  becomes: 
d (BCDm) 

- r_  _  -HCDmDE  +  (RCEDm* BCMSD)  (3.2) 

dt 
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and  rearranging  equation  E.2  produces 


d (ECDj) 

-  +fiCEjjDB 

dt 


FCLEm  = -  (3.3) 

FCMSD 

fchere  BCEEm  =  resident ia 1-comm ercial  distribution 

multiplier  for  fuel  m. 


RCt'.SD  =  market  sensitive  (free)  demand  for  the  BCD 
sector . 

BCDmDS  =  decline  rate  cf  BCD  demand  for  fuel  m. 
d (BCEm) 

- -  Ea+.e  Qf  change  of  BCD  demand  for  fuel  m. 

d  t 

The  application  of  equation  B.3  to  the  actual  data  will 
produce  the  desired  data  for  the  actual  distribution 
multiplier.  Therefore,  the  first  step  in  obtaining  the 
actual  distribution  multipliers  was  to  calculate  the  free 
demand  in  each  sector  from  the  collected  energy  demand  data. 
Although  no  problems  existed  with  calculating  the 
incremental  demand  portion  of  the  free  demand,  it  was 
necessary  to  use  the  earlier  assumptions  regarding  the 
consumer  commitment  liberation  rates  in  order  tc  calculate 
the  replacement  demand  portion  cf  the  free  demand .  Once  the 
free  demand  was  calculated,  the  decline  rate  cf  demand  for 
each  fuel  and  the  incremental  demand  for  each  fuel  were 
easily  determined. 

Tables  E.1  to  E*4  contain  the  actual  distribution 
nultipliers  as  calculated.  By  definition,  these  distribution 
multipliers  must  sum  to  1.0  for  any  given  time  period  and 
individually  they  can  never  be  less  than  0.0. 


- 

k#y  ! 

. 

. 

. 
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From  the  data,  several  parameter  estimations  were 
performed1  fcr  each  sector  assuming  different  equation 
structures.  Only  the  final  results  are  presented  here.  In 
the  instances  where  non-linear  structures  were  accepted  for 
the  distribution  multilpiier s  it  was  necessary  to  make 
adjustments  tc  the  equation  tc  insure  that  their  slope  is, 
at  every  point,  ncn-positi ve.  These  adjustments  were 
accomplished  ky  arbitrarily  cheesing  a  slightly  smaller 
price  ratio  than  that  which  yielded  a  minimum  fcr  the 
distribution  multiplier.  From  the  point  on  the  non-linear 
curve  corresponding  to  this  price  ratio  a  straight  line  was 
specified  which  had  an  equal  first  derivative  to  the  initial 
function  at  the  point  of  intersection.  Wherever  these 
adjustments  occured  they  will  be  presented  with  the 
parameter  estimation  results.  As  a  result,  the  equations 
described  below  for  the  distribution  multipliers  are  those 
which  were  used  in  the  'Canadian  Interfuel  Substitution 
Model 1 . 

RESIDENTIAI-CCF  KERCI AL  SECTOR 

ECIEWl  =  35.1764  -  .3437*FWB ATIC 
(2.98)  (.073) 

-  . 01176*TIME,  R2  =  .9123  (B.4) 

(.0015) 

if  EXE  AT  ICC. 42, 

ECECX1  =  -  17.9388  -  .  4 5 1  4* F XE ATIO  +  .5336*  (FX RATI  0) 2 
(2.42)  (.289)  (.  336) 

+  .G093*TIME,  R2  =  .7746  (3.5) 

(.0012) 


1  the  technique  used  was  ordinary  least  squares. 
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. 
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if  FXBATIO.  <12, 

RCEIX1  =  -18,0343  -  . C031 8*  (PXBATI0-. 42) 

+  .CG93*TIME  (B.6) 

The  best  estimate  for  the  equation  for  RCEDY1  was  found 

ty  t rial-and-er ror  and  is  equal  to  a  constant. 

ECDDY1  =  .44  ( B . 7) 

RCEDZ1  =  -19.3C65  -  .02  1  2* EZE ATIO 
(3. 384)  (.0248) 

+  .0101* TIME,  R2  =  .6509  (B.8) 

(.0018) 

In  every  instance  where  a  statistically  estimated 
equation  was  accepted,  all  of  the  estimated  parameters  were 
significantly  different  from  zero  in  a  statistical  sense 
(9011  level  of  confidence)  except  for  the  coefficient  of 
FZB ATIO  in  the  equation  for  RCDDZ1.  Despite  this  fact,  the 
coefficient  of  PZRATIQ  was  maintained  at  its  estimated  value 
for  use  in  the  model. 

INDUSTRIAL  SECTOR 

Since  the  coefficient  of  P ft RATIO  in  the  equation  for 
IHIEWl  was  net  significantly  different  from  zero  and  it 
possessed  the  wrong  sign;  it  was  ommitted  from  the  equation. 

IHDDW1  =  21.1355  -  .C107*!IME,  B2  =  ,5651  (B.Q) 

(3.863)  (.0020) 

if  EXRAT I0< . 40 , 

IHDCX1  =  -26.8776  -  . 6 5  1  9* F XR ATIO  +  .  7956*  (PXRATIO)  2 
(2.839)  (.348)  (.393) 

+  . 0 1 39*TIMF  ,  E2  =  .  8395  (B.10) 

(.0015) 


. 


. 


. 
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if  EXEAT TC>* 40, 

THE IX 1  =  -27.0111  -  . 0154*  (PXHATIO-.40) 

+  .0135* TIME  { B . 1 1 ) 

1  HEDY  1  =  -9.5911  -  . 01 2*E YRATIO 
(3.893)  (.061) 

+  . 005 G4* TIME  ,  R2  =  .1821  (B.12) 

(.  00  2) 

IHED21  =  12.5621  -  .  01  1 2*EZB ATIG 
(6.282)  (-044) 

-  .0064* TIME,  B2  =  .2061  (B.13) 

(.0033) 

Although  the  coefficients  of  PYEATIO  and  EZRATIO  in  the 
equation  for  IHEDY1  and  IHDDZ1  respectively  were  not 
significantly  different  from  zero,  they  were  net  removed 
frciB  the  equations  fer  use  in  the  model. 

TRANSPORTATION  SECTOR 

TREDW1  =  45.8318  -  . 271 8*EWRAII0 
(5.  537)  (.1  39) 

.  G232*1IME.  E 2  =  .7782  (B.14) 

(.  0  028) 

TREEY1  =  -47.6536  -  .  007 7*P YEATIO 
(5.799)  (.095) 

+  .024  8*TIME ,  E2  =  .  7395  (B.15) 

(.003) 

The  coefficient  of  PYRATIO  was  maintained  in  the  model 
at  its  estimated  value  despite  the  fact  that  it  is  not 
significantly  different  frem  zerc  at  a  90%  level  of 


confidence. 


. 


. 
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ELECTRICITY  GENERATING  SECTOR 
if  EW  E  AT IC< . 7  4 

ZEDW1  =  -16.157*1  -  4„ 9351*PWBATIO  +  3 . 29  3  1*  ( P  WE  ATI  C)  2 
(13.876)  (3.959)  (3.  032) 

+  .0095* TIME,  R2  =  .0764  (B.16) 

( .0073) 

if  EWR AT IO> , 74 

ZCDW1  =  -  18.04607  -  . 06 1 3*  (P WE A TIG- . 74 ) 

+  .CC93*TIKE  ( B . 1 7 ) 

if  E  XR  AT IG< . 66 

ZDDX1  =  1*1.2029  -  1  . 951 G*PXRATIO  +  1 . 4762  *  (PXFATI  C)  2 
(8.566)  (1.069)  (1.181) 

-  .  0  C  6  9* IIME ,  B2  =  .4129  (B.18) 

(.  0044) 

if  EXE  AT  10.66 

ZDDX1  =  1  3.5644  -  .  0024 1 6*  (PXPATIO- . 66) 

-  . 0069*11 ME  ( B . 1 9) 

if  E YR ATIO<<  1.21 

ZDDY  1  =  .5745  -  .  7 1  04  * F Y E AT 10 

(.299)  (.385) 

✓ 

+  .  2  92  2*  (P YE AT 10) 2  ,  E2  =  .4597  (B.20) 

(.123) 

if  EYEATIOI  .21 

ZCDY1  =  14.273  -  .0  032  8*  ( P Y RATIO-  1 . 2 1 )  (B.21) 

For  the  primary  demand  sectors  it  is  interesting  to 
note  that  only  the  distribution  multiplier  for  natural  gas 
responds  to  the  relative  fuel  price  in  a  non-linear  manner. 
This  has  been  interpreted  by  Khazzoom  to  mean  that  the 
ir.fluence  of  relative  price  changes  declines  as  the  relative 


price  increases. 
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TAELE  E  .  1 


CALCULATED  DISTRIBUTION  MULTIPLIERS  FOE  THE  E  ESI  DENT  IAL- 

COMMERCIAL  DEMAND  SECTOR 


NATURAL 

CRUDE 

YEAR 

CCAL 

GAS 

OIL 

ELECTRICITY 

T9  4  5 

0.4  7 

C  .  06 

0.32 

0.15 

1946 

0.41 

0,06 

0.36 

0.17 

1947 

0.46 

C  .  0  9 

0.15 

0.27 

1948 

0,42 

0.1 

0.22 

0,26 

1949 

0.34 

€.06 

0.3  6 

0.24 

1950 

0.3 

C  ,  07 

(5.39 

0.24 

1951 

C.  €8 

C .  05 

0.49 

0.37 

19  52 

0.26 

0.07 

0.38 

0.29 

195  3 

C.  18 

0.1 

0.45 

0.26 

1954 

0.23 

0,  1  1 

0,3  8 

0.27 

195  5 

0.  15 

0.1 

C  .  47 

0.28 

1  956 

0.09 

0,  19 

0,37 

0.35 

1957 

C  .  1 3 

0.2 

0.4 

0.28 

1958 

0.08 

0.2  3 

0.38 

0.3  2 

1959 

C  .07 

0,22 

0.3  9 

0.32 

1960 

0.09 

0,18 

0.36 

0. 37 

1961 

0,05 

0.2 

0.4 

0.35 

1962 

0.04 

0.1  5 

0.4  9 

0  „  3  2 

7  963 

0.03 

0,19 

0.45 

0.33 

1964 

0.0  1 

0.2  1 

0.48 

0.3 

1965 

0.01 

0.2 

0.32 

0,47 

1966 

0.01 

0.21 

0.3  7 

0.41 

1967 

0.0 

0.2 

0.34 

0.46 

1968 

0.01 

0.25 

C  .  33 

0.41 

1969 

— 

— 

- - 

1970 

— 

— 

- — 
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TABLE  B .2 

CALCULATED  DISTRIBUTION  MULTIPLIERS  FOR  THE  INDUSTRIAL 

DEMAND  SECTOR 


NATURAL 

CRUDE 

YEAR 

COAL 

GAS 

OIL 

ELECTRICII 

1945 

0.31 

0.04 

0.1  1 

0.54 

1 S46 

— 

- - 

- — 

— 

1947 

0.24 

C  .  04 

0.24 

0,48 

1948 

0.31 

0.03 

0.12 

0.53 

19  4  9 

0.32 

C  .  04 

0.09 

0.54 

1950 

0.23 

0.04 

0.  18 

0.55 

1951 

C  .  2  1 

C  .  03 

0.12 

0.6  4 

1952 

0.36 

0.03 

0.09 

0.52 

1953 

G.G8 

C.C8 

0.29 

0.55 

1954 

0,31 

C.  1 1 

0.  19 

0.39 

1955 

G .  25 

C.  1  1 

0.25 

0.39 

1956 

0.24 

0.12 

0.16 

0.47 

1957 

0.16 

C .  1  3 

0.15 

0.56 

1958 

0.  13 

0.2  7 

0.2 

0.4 

1959 

C.  15 

0.17 

0.17 

0.52 

I960 

0.17 

0.27 

0.29 

0.26 

1961 

0.17 

0.2 

0.24 

0.39 

1962 

0.18 

0.2 

0.28 

0.34 

196  3 

0.18 

0.18 

0.15 

0.49 

1964 

0.17 

0.2 

0.  19 

0.44 

19  6  5 

0.08 

0.25 

0.23 

0.4  3 

1966 

0.12 

0.24 

0.26 

0.38 

1967 

0.05 

0.2  6 

C.34 

0.35 

1968 

0.02 

0.29 

0.17 

0.53 

1969 

— 

— 

— 

— 

1970 

— 

— 

- - 

— 
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CALCULATED 


LADLE  B.3 


DISTRIEUTICK  MULTIPLIERS  FOE  THE  TRANSPORTATION 
DEMAND  SECTOR 


CRUDE 


YEAR 


1945 
1  946 

1947 

1948 

1949 

1950 

1951 

1952 
19  53 
1  954 

1955 

1956 

1957 

195  8 
1859 
1  860 
1961 
1  962 
1  963 
1  964 
1965 

196  6 
1867 
19  6  8 
1  96  9 
1  970 


CCA  I 


0.52 
0.59 
0.38 
0.46 
C  .2 
0.41 
0.5 
0.46 
0.35 
0.33 
0.22 
0.04 
C  .04 
0.01 
G.C3 
0.02 
C.03 
0.02 
0.05 
0.04 
0 .01 
0. 

0  . 

0. 

0 .01 


CIL 


0.48 
0.41 
0.6  2 
0.54 
0.8 
0.5  9 
0.5 
0.54 
0.65 
0.67 
0.7  8 
0.96 
0.96 
0.99 
0.97 
0.98 
0.97 
0.98 
0  .95 
0.96 
0.99 
1.0 
1  .0 
1.0 
0.99 
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T  A  E  LE  E .  4 


CALCULATED  DISTRIBUTION  MULTIPLIERS  FOR  THE  ELECTRICITY 

GENERATING  SECTOR 


NATURAL  CRUDE 


YEAR 

CCAL 

GAS 

OIL 

1945 

0.73 

0.1 

0.17 

1946 

0.7  5 

C  .09 

0.16 

1947 

0.74 

0.  1 

0.16 

19  48 

0.79 

C  .07 

0.14 

1949 

0.67 

0.  1  5 

0.18 

1950 

0  .  C  4 

0.54 

0.42 

1951 

0.54 

G  .  2  2 

0.24 

1952 

0.6 

C  .2  1 

0.19 

1953 

0.35 

0.38 

0.26 

1954 

-- - 

— 

— 

1955 

0.61 

0.2 

0.19 

1956 

0.57 

C  .  3  2 

0„1  1 

1957 

— - 

— 

— 

1956 

- - 

— 

1959 

— 

— 

1960 

0.42 

0.42 

0.16 

1961 

0.18 

0.  1  1 

0.08 

1962 

0.64 

C.  1  3 

0.02 

1963 

0.77 

0.07 

0.16 

1964 

0.62 

0.2  4 

0.14 

1965 

0.48 

0.32 

0.2 

1966 

0.63 

0.11 

0.25 

1967 

0.56 

C.  1  4 

0.3 

1968 

■ - - 

— 

—  — 

1969 

-  —  - 

— 

1970 

—  — 

— 

— 

APPENDIX  C 


ECDEI  LISTING 

This  appendix  contains  a  listing  of  the  "Canadian 
Interfuel  Substitution  Model".  The  computer  language 
employed  is  hr.cwn  as  BYNAMC  (21). 

The  exogenous  variables  in  the  model  as  it  is  listed 
here  are  those  that  were  assumed  for  CASE  III  of  the 
projections. 
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